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|. INTRODUCTION 


A guided missile is one which is usually fired in a direction 
approximately towards the target and subsequently receives steering 
commands from the guidance system to improve its accuracy. 

In recent years, the application of Bank-to-Turn (BTT) guidance to 
tactical missiles has generated conciderable interest, motivated by certain 
advantages that a Bank-to-Turn control configuration can offer. 

First, against high performance threats, there is a need for defensive 
missiles to develope increasingly higher lift accelerations, which requires 
banking maneuver to properly direct the control vector. 

Second, the need for greatly increased stand off ranges has led the 
design of tactical missiles towards airbreathing propulsion systems such as 
the ramjet. These designs generally have configuration geometries that are 
not cruciform, due to exposed inlets beneath the vehicle. As a result, there 
are often stringent limits on the sideslip angle that can be developed during 
engine operation, and this too, dictates some types of Bank-to-Turn control 
scheme. 

This thesis investigates the controllability of a Bank-to-Turn missile, 
when one or more control surfaces are damaged, at fixed flight conditions 
(i.e. constant altitude, Mach number and missile weight and inertias). 

The approach is to design a control reconfiguration module to be placed 
between the flight control laws and the control surfaces. In the face of 


missile damage the reconfiguration module would reallocate control 
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commands to the unimpaired contro} surfaces so that flight critical pitch, 
yaw, and rolling moments would be preserved to the maximum extent 
possible. Using the reconfiguration module control laws, the existing flight 
control laws would not be altered. 

In the first part of the thesis the control surfaces conventions, 
definitions and the manoeuvre by banking are reviewed. 

The second and third part reviews the design of the longitudinal and 
lateral channels and involves the uncoupled pitch, roll and yaw channels for 
the circular and elliptical airframes of a B.T.T. missile. The resultant 
continuous open loop designs are analyzed and found to be in accordance 
with the desired specified requirements. Applying the convention contro! 
theory, the corresponding modified longitudinal and lateral autopilots are 
obtained, analyzed and compared with the original system. The effect of 
damage to one or more control surfaces, on the controllability of each 
system, is obtained, analyzed and compared with the unimpaired models 
respectively. 

Finally, the combined linear open loop pitch and roll channel autopilots, 
for the circular airframe configuration, are investigated. The effect of 
damage to one or more actuators, on the control of the combined unimpaired 
system is obtained and analyzed. The most severe type of impairment 
considered in this thesis, i.e., the hard over type, is extensively investigated 
in Chapter Y. Corrective actions are then developed and presented. 

The robustness.of the complete system is analyzed in chapter five. 
Conclusions and recommendations for future study are then given. 
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The investigation and analysis in all the above cases was performed 
using the existing at Naval Postgraduate School (NPS), Optimal Systems 
Control Program (OPTSYS) for the continuous system, Robustness Design 
Program (POPLAR) for the robustness of the system, and EASYPLOT 
program for the plotting of different kind of impairments and comparisons. 
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ll. GUIDED MISSILE CONTROL SURFACES CONVENTIONS AND MANDEUYRE 


In this chapter, before going into mathematical details concerning the 
motion of a guided missile in space as a result of guidance commands, some 


definitions and discussion concerning the control surfaces are given. 


A. GENERAL 
The task of a Missile Control system is to manoeuvere the missile 
quickly and efficiently as a result of the following two control singnals: 


(a) Detection whether the missile is flying too high or too low, or too 
much to the left or right. 


(b) Measure the deviations or errors and develope commands to the 
control system to reduce these errors to 2810. 


B. CONTROL SURFACES DEFINITIONS 

Suppose that the control surfaces 1 and 3, as is shown in the Figure 
2.1, are mechanichally linked together such that a servo must impart the 
same rotation to both. Then these surfaces are essentially elevators. 
Positive elevator deflection produces a negative force in the 
2-direction and an anti-clockwise moment about the y-axis. 

If now surfaces 1 and 3 each have their own servo it is possible for 


them to act as ailerons. Positive aileron deflection produces an anti- 





clockwise moment about the x-axis. 
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If looking in the direction y one surface is rotated 5° clockwise and 
the other surface 5° anti-clockwise then a pure couple is imparted to the 
missile about the fore and aft axis and this will tend to make the missile 
roll. Such control surfaces are now called ailerons, and by doing the same 
thing to control surfaces 2 and 4, the turning effect of the ailerons can be 
doubled. 





Figure 2.1 Control surfaces looking from rear of the Missile. 


If the aerodynamics are linear i.e., the normal forces are proportional 
to incidence then the principle of superposition applies. Commmands for 
elevator rudder and aileron movements can be added electrically resulting in 
unequal movements to opposite control surfaces. Positive rudder 
deflection produces aé positive force in the y-direction and negative 
moment about the z-axis. In this way we have the means to control roll 
motion as well as the up-down i.e., pitch motion and left-right i.e., yaw 


motion. 
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С. CONTROL SURFACE CONVETION 
As is shown in the Figure 2.1 deflections 5,, 55, 5z, 54 are defined as 
positive if the leading edge of the control surface is rotated clockwise as 
one looks outwards along the individual hinge axis. The total aileron 
deflection is defined as: 
Aileron deflection 5=1/4 (5, +55+52+5,) (II.C- 1) 
or (+5+ן2)5/ן‎ or 1/2(825*6,4) (П.С-2) 


if only two surfaces act differentially. 


D. BANK-TO-TURN CONTROL SURFACES CONFIGURATION 
The same kind of analysis can be applied for the Bank-to-Turn missile 
configuration, shown in Figure 2.2. Note that the arrow on the control 


surface indicates the direction of positive leading edge deflection. 


Tail Tail Incidence 


omm 51 5৫ [5৩54 


+ 


ממממוט 


Figure 2.2 B.T.T Aerodynamic Sign Convention. 





The maneuver plane is the xg-zg plane. The downwards movement is 


defined as positive. This is the configuration investigated in this thesis. 
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The B.T.T. missile configuration is tail-controlled using four identical 
control surfaces. The control surfaces are located flush with the body base 
with a +30 degrees dihedral. With tail control configurations the pitching 
moment or control effectiveness is quite high in view of the fact that the 
tail moment arm is generally fairly large. Also there is no adverse 
downwash from deflection of the forward surfaces, in comparison with the 
canard control configuration. A more complete model of the B.T.T. circular 
confinguration is also shown in Figures 3.1 and 3.4, taken from Arrow, А. 
[Ref. 1]. 


E. MANEUVER BY BANKING. 

The Conventional method of altering course to the right or left in a 
missile or aircraft is to use the ailerons to bank i.e, roll by an angie 9, as is 
shown in Figure 2.3. If the lift force is increased slightly by the use of the 
elevators so that the vertical component of lift equals the weight then 
there is a horizontal component of lift equal to the total lift times sin 9. It 
is this component of lift which causes the flight path to change. An exact 
analysis of the maneuver is not simple as the airflow over the outer control 
surface is faster than that over the inner control surface due to the rate of 
turn. In practice a small amount of rudder is applied in an attempt to make 
the general airflow directly along the fore and aft axis of the cordinated 
B.T.T. missile and in the plane of the control surfaces. In this condition 
there is no 'sideslip and hence no net sideforce. This is the preferred 
method of manoeuvring since lifting forces are most efficiently generated 


Perpendicular to the control surfaces; the lift-to-drag ratio is a maximum 
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in this condition. Even so there is some additional induced drag when 
manoeuvring. An additional manoeuvring force is obtained by increasing the 
incidence and this means that the normal force is inclined backwards (0 the 
velocity vector by the angle of incidence. There is a useful component equal 
to the normal force times the cosine of the angle of incidence and a 
component proportional to the sine of the incidence opposing the velocity 
vector. This “induced drag’ will always occur whenever a lifting force is 


produced irrespective of the method of manoeuvre. 


vertical 
component 
of lift 







horizontal 
component 
af lift 

weighit 


centrifugal 
force 





Figure 2.3 Forces Acting on a Banking Missile and Coordinated Missile 
Motion for coordinated B.T.T. Control Policies. 


Because it is desired that the B.T.T. missile avoid negative angles ої 
attack, the missile is forced to roll about its velocity vector when the 
desired manouvre direction is in the negative angle-of-attack direction. 
Although some negative angles of attack may be possible with certain 


missiles. 
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Coordinated motion or zero sideslip angle is achieved by directing the 
body fixed pitch axis of the missile at the missile velocity vector so that 
there is no component of missile velocity vector along the body fixed yaw 
axis of the missile. The Figure 2.3 (b) shows the attitude of the missile 
body with respect to its velocity vector ¥. When an upward command is 
applied, i.e., ф = 0, the missile body moves upward with its pitch axis 
directed at the velocity vector until it reaches manoeuver level or angle of 
attack. No roll motion 1s required to maintain coordination of this 
manoeuver. The circular airframe configuration of this thesis is a planar 
configuration which is symmetric about its wing planes and therefore have 


two perfect directions which are normal to the planes of its wing. 
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LONGITUDINAL MODEL METHODOLOGY‏ .וו 


A. INTRODUCTION 

in this Chapter the design of the continuous uncoupled pitch channel 
(longitudinal) autopilot is reviewed, as was analyzed by Arrow, A. and 
Karaiskos, |. [Ref. 1 and 2] A modified longitudinal autopilot is then 
designed with the four control surfaces split into independent components, 
for the circular configuration model shown in Figure 3.1. 

Then a comparison between the two systems (original and modified) 
shows that are identical. The transient response plots and their analysis 
demonstrate the performance of the modified pitch channel. 

Having established the validity of the modified system one then uses 
the multiple-input multiple-output control theory and state-space 
representation [Ref. 3), to investigate the performance of the modified pitch 
channel autopilot for various classes of damages. 

For this thesis, two classes of impairments are considered namely, 
locked or hard over and inoperative with the control surfaces centered. The 
plots of a modified pitch autopilot without impairment are used as a basis 
for comparison, the impairment induced effects on normal acceleration, rate 
and tail incidence are assessed. 

À general block diagram, of a B.T.T. autopilot is shown in Figure 3.2 
taken from [Ref. 1]. Inertial acceleration commands are applied in polar 
coordinates (i.e, magnitude of the command ( mc ) applied to the pitch 
autopilot and the direction ( 9c ) is applied to the roll autopilot. Achieved 
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manoeuver plane or inertial acceleration in rectangular coordinates (i.e, Tz 
and jy) is determined by resolving achieved body-fixed accelerations ( i.e., 
Nz and (y) through missile roll angle (99). 

Before going into the details concerning the various kinds of damages 
it is prefered to make here clear some terms used in this work. 

ORIGINAL SYSTEM is called the unmodified pitch channel 
autopilot presented in [Ref. 1]. 

UNIMPAIRED SYSTEM is called the modified, i.e., with the split 
control surfaces system, and without any kind of damage. 

IMPAIRMENT, HARD OVER OR INOPERATIVE are terms expressing 


a specific kind of damage to one or more control surfaces. 


B. UNCOUPLED LINEAR PITCH CHANNEL (ORIGINAL) 
1. Aerodynamic Model 
The linear pitch channel dynamic model which is used for the 


circular airframe, is shown in Figure 3.3, with the following requirements : 


а. In actuator command branch, the high attenuation frequency must be 2 
15 db at 100 rad/sec and zero angle-of-attack. This requirement 
limits the autopilot speeds of time response. 


b. Relative stability with gain margins 26 db, phase margins 2 0 
degrees with a goal at 12 db and 50 degrees. 


c. An acceleration time responce with 63% time constant of 0.5 seconds 
for a step command of acceleration at the flight condition of interest, 
overshoot < 10%, and Zero steady-state error to reduce variations of 
guidance navigation gain. 
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The linear pitch channel dynamic model which is used, is shown in 

Figure 3.3. This Figure is taken from [Ref. 1]. The model is linearized for 
Stability studies and the following three assumptions were made : 

a. Manoeuver plane is the plane xg - 2g , as is shown in Figure 3.4 in 


which the arrow indicates the direction of positive leading edge 
deflection. 


b. The missile is trimmed in pitch, at fixed values of angle of attack (8), 
pitch angular rate (q), and pitch tail incidence (Sp), that is M,= 0. 


c. The roll rate (P) of the missile is constant. 

Using the foregoing assumptions, and the approximation that the 
variation in forward velocity (u), is equal to zero a simplified set of 
longitudinal equations of motion are constructed applicable to the short 
period of motion. These equations describe the two degree of freedom, short 
period longitudinal mode. 

Linearized aerodynamic derivatives for M = 3.95 and Altitude = 60 
KFT are provided in Table 3.1, as is given in (Ref. 1]. 

2. Pitch Control Law 

The longitudinal control law model for the circular airframe was 
developed from the block diagrams and information provided in [Ref. 1], and 
is shown in Figure 5.5. Lead-lag filters were used in the system to shape 
the output response, provide adequate gain and phase margins, and finally to 
prevent guidance noise saturation problems. The guidance requirement of 
zero steady-state error is satisfied by using an integrator іп the 
acceleration error branch of the circulare control law. The rate error 
conpensation determines the high frequency attenuation and was used to 


minimize the effect of aerodynamic variations on acceleration time 
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response. The acceleration error compensation determines the acceleration 
time response. 

The blocks in the pitch control law of Figure 3.5 for the circular 
model represent the basic transfer functions in which a normal acceleration 
command (nz,) equal to 1 gee, is applied to the pitch control low which uses 
measurements of: 

a. Missile body pitch angular rate (q) and 
b. Pitch normal acceleration (ny) 


to determine the required actuator comand (6рг). Тһе actuator is modeled as 


a first order lag at 188.4 rad/sec and is shown in Figure 3.2. 
5. Basic Transfer functions of the Aerodynamic Model 
The equations which are represented by the block diagram of the 
aerodynamic model, of the uncoupled pitch channel are presented below and 
taken from [Ref.2] . 


a. The pitch aerodynamic transfer function for the pitch 
angular rate (q) about the Yp axis: 


q -1.336 1s- 0.15 
— z= ——— Ñ  —c eG) (111.B.3- 1) 
Sp 22.545 х 1075 s2 + 3.5634 א‎ 1055 + 1 


b. The pitch aerodynamic transfer function for the achieved 
maneuver normal acceleration (nz) in the Zg direction: 


n, 106.074 x 10-3 s2 - 18.82 


— -— T — gra) — (HIB3-2) 
б 22.545 х 10-3 52 + 3.3634х 10-55 + 1 
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Figure 3.4 Aerodynamic Sign Convention and Axes. 
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TABLE 5.1 


LINEARIZED AERODYNAMIC DERIVATIVES (M=3.95) 
FOR CIRCULAR AIRFRAME 


8509 8510 8520০ 
-0.065 -0.082 -0.111 
201025 -0.019 -0.003 
0 -0.009 -0.020 
0.021 0.022 0.028 
-0.050 -0.053 -0.062 
0 -0.016 -0.038 
0 -0.009 -0.022 
0 0.018 0.044 
0.051 0.035 0.044 
9.15 0.17 0.22 
0.04 0.04 0.05 
-0.060 -0.065 -0.118 
-0.080 -0.095 ₪ - 


58 


Reference Center of Gravity (C.G) at 0.6 


4 Analysis for Circular Airframe 


The analysis of the uncoupled pitch channel autopilot is based on 
the time responses and Bode plots of maneuver plane acceleration, body 
angular rates and tail incidence angles, obtained by OPTSYS control program 
of Naval Postgraduate School (NPS), which requires the differential 
equations to be in the state-space variable format: 


Х=ЕХ +60 (III.B.4-1) to 


define the system. In this system the F matrix is the open loop dynamics 
matrix (system or plant) and the G matrix is the control matrix. The 
variable assignements are X as the state vector and U as the control input 
vector. It follows that X is the time derivative of X. 

Using state-space representation of a system, in which the 
forcing function involves derivatives, the following equations are derived 
and can be modeled as a tenth order system, for the uncoupled pitch 


channel : 
Х1= -150 Х, +150 Хь + 705.25 65, (111.8.4-4) 
X, = Xz - 1.3531 ₪ + 1.3531 * (111.B.4-5) 
Xs = -5 Xs + 6.572 na, - 6.5?2 X, (111.8.4-6) 
Х4 = Хз - 59.2637 6, (111.B.4-7) 
Xs = -149.186 х10-3 Xs -44.355? X4 + 1.7887 8, )!!1.8.4-6( 
Xg z X? - 701.92х10:3 6, (111.B.4-9) 
Хә = -149.186x10-3 X; - 44.3557? X& - 1043.36 6, — (I11.B.4-10) 
Š, = 7188.4 6  188.4 бр, (111.8.4-11) 
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áp, = -0.143 dp, -2.2308 Х, % 2.2308 Х, - 540.58х10:3 Х; 

* 460.839x1073 na, - 460.839х10-3 Х, + 340.58 

x1073 Xs - 20.184 6, (111.B.4- 12) 
fy = -( Xy - 887.124 8, * 886.422 бр ) (111.B.4-13) 


where the state variables are: 
e X, X? .pitch angular rate state variables 


€ X; X; : achived body acceleration state variables 
ex . output of acceleration filter 
е бр: : input command in actuator network 
өб . tail incidence angle 
Expressing the above equations ав а set of first order matrix 
in the format of (III.B.4-1) the plant and control matrices are: 


41-150 0 0 0 0 0 0 0 0 150 
1.353 0 1 0 0 0 0 0 0 0 
-6.572 3 0 0 0 0 0 оо 

0 о 0 0 l 0 О -59.26 0 0 

F= 0 0 0 -4456 -0.15 0 0 1.789 0 0 
0 о 0 0 0 0 | 9 0 0 

0 0 0 0 O -44.36 -0.1492 -1043 0 0 

0 0 0 0 0 0 0 -188.4 188.4 0 
-0.461-2.251-0.5406 2.251 0.5406 0 0 -20.18 -0.143 0 

Ó Ó Ó 0 0 0 -1 887.1 -886.4 0 


(III.B.4- 14) 


The control matrix 6 is following: 


-1.3531 

6.572 
0 

6= 0 (III.B.4- 15) 

0 
0 
0 

0.4608 
0 
where: U = nz, as control input to the system. (III.B.4- 16) 


According to the [Ref. 1], using a step input function which 
represents the “1 gee Command” and a trim angle of attack a,=0, the time 
response plots of the state-variables of interest were obtained from the 
OPTSYS program and are represented at the Figures 3.6, 3.7 and 3.8. Typical 
performance criteria used to characterize the transient response to the one 
gee unit step input include overshoot, delay time; rise time, and settling 
time. The above-mentioned criteria are defined with respect to the step 


response: 


e /eximum Overshoot. The maximum overshoot is defined as the largest 
deviation of the output over the step input during the transient state. 
The amount of maximum overshoot is also used as a measure of the 
relative stability of the system. The maximum overshoot is 
represented in this thesis as a percentage of the final value of the 
step response: that is, 


max. overshoot 


percent maximum overshoot = —————————- x 1003 (IIL.B.4-17) 
final value 
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e )Dejay time The delay time is defined as the time required for the step 
response to reach 508 of its final value. 


е se time The rise time is defined as the time required for the step 
response to rise from 10% to 903 of its final value. Frequently the 
rise time is also refered to as the time constant of the system. 


e Settling time The settling time is defined as the time required for the 
step response to decrease and stay within a specified percentage of 
its final value. For this thesis is specified arbitrarily at 0.012. 


The four quantities defined above give a direct measure of the 
transient characteristics of the system to the step input: 
a. The pitch normal acceleration (nz ro Nz), Figure 3.6, g's vs. 
time which has: 
(1) A05 rise time or time constant. 
(2) А тіпог 0.005 steady state error and 2.2% overshoot. 
(3) A 0.28 sec delay time and 2 sec settling time. 
b. The pitch angular rate (q or Q), Figure 3.7, rad/sec vs. time, 
which has: | 


(1) The maximum body angular rate (-0.17 rad/sec) occurs at 0.28 second 
when the tail incidence is at 0.02 rad angle of deflection. 


(2) From the moment at 0.28 sec the rate increases and reaches at 2 sec 
its steady-state value (-0.01rad/sec). 


(3) The -0.01 rad/sec steady-state constant rate is due to the fact that 
one has a type zero system and therefore has a finite error to a step 
input. 


с. The pitch tail incidence (êp or DP), Figure 3.8, rad vs. time, 


which has: 


(1) Actuator deflection angle 0.05 rad. 
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(2) Initially and within the 0.28 sec delay time a small transient occurs 
with an insignificant fluctuation. 


(3) An overshoot 1.3% and a 2.0 sec settling time. 
The above analysis has revealed that the perfomance of the system is 
according to the previous mentioned requirements. 
Also should be mentioned here that according to the Eigenvalues found 
executing the OPTSYS program, the continuous open loop system is stable, 


since the s-plane poles are: 


Sı = -159.724 + ] 19.9921 


52 = -159.724 - |18.9921 
S3 = -8.29048 + 2 
S4 = -8.29048 - |8.05932 
55 = -3.75925 + 353 
56 = -3.75925 - |2.51463 


Sp = -0.14393 + jo 
Sg - -0.00000425 % 10 
Sg = -0.04593 + j6.65959 
510: -0.04593 - |6.65959 
5. Design of the Modified Uncoupied Pitch Autopilot 
To study reconfigurable flight controls the missile model should 
be capable of using the set of four available control surfaces to produce the 
required forces and moments. 
To achieve this each control surface ( 84, 82, 83, 64 ) is considered 


as an independent component. 


The equations are then coupled so that a component of control 
surfaces used as a single or pair, will produce the appropriate moments. For 
example ô; deflected at an certain angle, while the other three remained at 
a zero angle, will produce primarily a rolling moment, and deflected 
together with 52 no roll effect is present, as we will see in Chapter V. 

Taking under consideration the conyention control surfaces theory 
of Chapter II., for the B.T.T. missiles, and for the configuration shown in the 
Figure 3.4, one can easily show, for the continuous pitch autopilot that the 
input command in the actuator network (65,) can be split up into independent 
components, using the following equation which compute the deflection 
inputs to Equation ( III.B.4-12 ). (Also refer to Figures 2.2 and 3.4 which 


show the control surface positions and corresponding positive deflections). 
+Špc = (-64 +2 -63 +64 ( 2 )111.8.5-1( 


Replacing the input vector in Equations (III.B.4- 11), (III.B.4-12) 
and (III.B.4-13) with the above Equation and applying the superposition 


principle, gives the following modified airframe equations: 


I 
бр = -188.46p * 188.4 — (- 8, * 82 - 83 * 64) )111.8.5-2( 
2 
I | 
—N- $; + $, < $; + 64) = -0.143 -- (-61 + 82 7 63+ 64( B 2.2308Х> 
2 2 
+ 2.2308Х4 - 340.58x10-3 X; + 460.839x10-3n,, 


- 460.839х10-5Х, % 340.58х10-5Х5 - 20.18465р 
(III.B.5-3) 


| 

fo, = - (X; - 8B7.124)5p + 886.422—{-8, * 82 - 53 + 84) (Ш.В.5-4) 

2 

In Figure 3.6 the modification of the commanded tail incidence is 
presented. 










ACTUATGR 


ACTUATOR р 


(ם) 


Figure 3.6 (a) Original System; (b) Modified System. 


Rearranging the above equations and using the transfer functions 
of the system for the uncoupled pitch channel, Equations (III.B.4-4 to 
III.B.4- 13), the following equations are derived and can be modeled in a 
thirteenth order system for the same channel: 


X, - -150X, + 1501; )111.8.5-5( 
גא‎ = Хз - 1.3531п,, + 1.3531Х, )111.8.5-6( 
Хз = -5Х:; + 6.572п,. - 6.572% )111.8.5-7( 


%4 = 5 - )111.8.5-8( 


Х5 = -149.186х10-3Х5 - 44.3557Х4 + 1.78875р (III.B.5-9) 
১6 5 ১7 -701.92%10-38 (Ш.В.5- 10) 
X; < -149.186х10-35Х; - 44.3557Хұ- 1043.365p )111.8.5-11( 
Sp = -188.48p - 94.25, + 94.255 - 94.253 + 94.25, )111.8.5-12( 


$, 2 -0.1438, * 1.1154X2 - 1.1154X4 * 0.17029X; 

0.23042X, - 0.17029%8 + 10.0296 (Ш.В.5-15)‏ + 2ח0.25042- 
1.1154Х; % 1.1154Х, - 0.17029Х;‏ - ;0.1435- - 855 

* 0.23042na, - 0.23042X, * 0.17029Xs - 10.0925 (Ш.В.5-14) 
z = -0.14384 * 1.1154X2 - 1.1154X4 % 0.17029Х; 

- 0.23042na, * 0.23042X, - 0.17029Xs * 10.0298p — (III.B.5-15) 
$4 2 -0.14384 - 1.1154X4 * 1.1154X4 - 0.17029X; 

* 0.23042nz, - 0.23042X, + 0.17029Xs5 - 10.0298 (Ш.В.5-16) 
r> = -X, * 887.1246p * 443.2118, -443.2118, * 443.21 18; 


- 443.2118, (III.B.5- 17) 


We thus have a multiple-input-multiple-output (MIMO) system 
with the following order: 
е: 13xi, the time derivative of X matrix. 
eF : 13x13, the open loop dynamics matrix or plant. 
@ X : 13x1, the variables assignements matrix. 
e 6 : 13x5, the control distribution matrix. 
el : 5x1, the control input matrix. 
The above plant control and output matrices are presented at the 
Appendix A of this thesis. 


6. Analysis and Comparison with the Original System 

Applying the same step input function as in section III.B.4, which 
represents the “1 gee command’, and zero degrees deflection command for 
the four control surfaces (i.e., 645 655 635, 64, - 0) the time response plots 
of the following state-variables of interest are obtained and presented in 
the Figures 3.10 through 3.13. 

The performance of the system ( N; q, б; ) as shown in Figures 
3.10 to 3.12 is the same as that shown in Figures 3.6 through 3.8. 

In Figure 5.15 is shown the time response of the control surfaces 
deflection angle, for the input of an acceleration equal to 1 gee command. 

Deflections ( 84, 82, 63, 64 or D's ), Figure 3.13, rad vs. time, 


which has: 


(1) The 5, and 85 have the same sign (negative), as it was expected 
according to the sign convention of Figure 2.2. 


(2) The 62 and 84, have the same response as above but with a positive 
sign. 

(3) Insignificant overshoot. 

(4) А 2 seconds settling time. 

In the same Figure 3.13 the deflection of the control surfaces is 
illustrated. The presented deflection is the required position of each 
surface in order to achive the commanded one gee acceleration (downwards). 

We procceed now to the analysis of the impaired control surfaces 


for the modified system. 
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Figure 3.7 Pitch Normal Acceleration vs. Time; Uncoupled Pitch Channel; 
Continuous Open Loop System; a, = 0; Step Input = 1 45. 
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3.8 Pitch Angular Rate vs.Time; Uncoupled Pitch Channel; 
Continuous Open Loop System; a, =0; Step Input - 1045. 
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Figure 3.9 Pitch Tail Incidence vs. Time; Uncoupled Pitch Channel; 


Continuous Open Loop System; a,=0; Step Input = 1 05. 
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Figure 3.10 Pitch Normal Acceleration vs. Time; Uncoupled Pitch Channel; 
Unimpaired Open Loop System; 4,=0; Step Input: n,,=1! g's and 
615 625 63 546 - degrees. 
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Unimpaired Open Loop System; a,=0; Step Input: n,,=1 g's and 
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Figure 3.13 Control Surfaces Deflections vs. Time; Uncoupled Pitch Chennel; 
Unimpaired Open Loop System; a,=0; Step Input: n,.=1 g's and 
610 bony ӛзе бас - 0 degrees. 
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C. CONTROL EFFECTIVENESS DUE TO THE CONTROL SURFACES 
IMPAIRMENT AND COMPARISON 


The time response plots of the system are obtained in this section 
for two types of damages. The first type is the locked or hard over 
surface at a small angle and the second type includes the inoperativity of 
one, two or three actuators. The second type of damage could also 
represent the destruction of the control surface(s), i.e., shot away. 

Note that the first phase of the design approach for the CBTT 
autopilot was to design each channel independentiy with all coupling 
between channels removed. Sufficient high frequency attennuation was 
added for actuator and missile design so that the resulting missile body 
angular rates and control surfaces motion would represent a practical 
missile design. A relationship was established among the relative speeds 
of response of the uncoupled channels in order to meet CBTT requirements. 
The cross-coupling due to control surfaces impairment between the 
longitudinal and lateral motion of the missile is ignored when each control 
channel is treated independently. Cross-coupling of the controls is thus 
considered in Chapter V. 

1. Hard Over One Control Surface 

One control surface (5,), was selected to be deflected or hard 
over, i.e., it remained at a fixed angle while the other three control 
surfaces have a variable deflection angle. 

The control input was a step input of one gee in the normal 
acceleration as in the original and unimpaired models. For the damaged 


control surface (6) a hard over deflection of 0.2 rad was considered. The 
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time response plots.of the state-variables of interest were obtained and 
analyzed with the performance criteria as defined in section (III.B.4): 
a. Pitch normal acceleration (n, or Nz), Figure 3.14, g's vs. 


time, which has: 


(1) Initially and with a fast response, the missile body tends to go 
upwards reaching a 0.6 gee acceleration. This is caused by the fact 
that at the beginning of the command the damaged surface was 
deflected at a +0.2 rad angle of deflection while the other three 
were at their steady position. Then the missile goes downwards 
until the -1.25 gee acceleration is achieved. This is due to the 
deflection of the unimpaired surfaces which are commaned in such a 
way as to counteract the damaged (5,) surface. From this point (0.2 
sec) the acceleration rises constantly in order to achieve the one 
gee command. 


(2) Arise time 0.75 seconds and a negligible overshoot. 
(3) A delay time about 0.5 sec 
(4) Steady-state error 0.005% which is insignificant. 
b. The pitch angular rate (q or Q), Figure 3.15, rad/sec vs. time, 
which has: 


(1) Initially and due to the fact that ô, is locked at a positive angle of 
0.2 rad, to achieve the acceleration response, the required body 
angular rate the first tenth of the second increases very fast to 
0.48 rad/sec. Then it begins to decrease crossing the zero rate at 
0.2 sec, corresponding to the max decceleration (-1.25 gee's). Then 
it changes direction for the next 0.2 sec until the value of -0.58 
rad/sec is achieved. From this point the rate is increased again 
until the constant value of -0.01 rad/sec is achieved at the 
settling time 1.4 sec. 


(2) The initial increase of the rate is due to the produced pitch moment 


caused by the deflection of the damaged control surface while the 
other three are at their steady position. 
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c. Pitch tail incidence (8, or DP), Figure 3.16, rad vs. time 


which has: 


(1) To achieve the acceleration response and due to the fact that 8, is 
locked at an angle of 0.2 rad, initally more control surface 
deflection is required ( -0.13 rad). At 0.2 sec, when the 
acceleration begins to increase and the angular rate is zero, the tail 
incidence has the positive angle of 0.05 rad and a fluctuation begins 
before being stabilized at a constant angle of 0.06 rad. 


(2) Settling time equal to 1.8 sec. 
d. 08161800075 (519 ê>, êz, 64 or D's ), Figure 3.17, rad vs. time. 
This Figure illustrates the time response of the deflection 

of the four control surfaces,when one of them (54) is locked at an angle of 
0.2 rad, in order to achieve the acceleration response. Initially while 62, 
65 and 54 control surfaces have zero angle of deflection, 84 is already at 
the 0.2 angle. Then 62 and $4 are deflected at a positive angle 0.1? rad and 
63 at a negative 0.1? rad, according to the sign convention shown in Figure 
2:2. | 

À comparison between Figure 3.1? and 3.15 of the unimpaired 
system, shows that to achieve the acceleration response the impaired 
system requires more control surface deflection. 

Finally to cover completely this class of impairment, the case of 
a negative angle (-0.2 rad) at the Š, surface is examined. The time 
response of the same state variables are obtained and analyzed. 

а. Pitch normal acceleration (п, ог №2), Figure 3.18, g’s 95. 
time, which has: 


(1) Initially and with a fast response, the missile body tends to go 
downwards reaching a 0.6 gee acceleration. This is caused by the 
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fact that at the beginning of the command the damaged surface was 
deflected at a -0.2 rad angle of deflection while the other three 
were at their steady position. Then the missile goes upwards until 
the +1.75 gee acceleration is achieved. The produced overshoot is 
due to the deflection of the unimpaired surfaces which are 
commanded in such a way as to counteract the damaged (5,) surface. 
From this point (0.3 sec) the acceleration deacreases constantly in 
order to achieve the one gee command. 


(2) A small 0.18 second rise time (reason for the big overshoot). 
(3) A delay time about 0.1 sec 


(4) Steady-state error 0.005%. 

b. The pitch angular rate (q or Q), Figure 3.19, rad/sec ys. time, 
which has: 

(1) Initially and due to the fact that 5, is locked at a negative angle of 
-0.2 rad the required body angular rate the first tenth of the second 
decreases very fast to 0.68 rad/sec. Then it begins to increase 
crossing the zero rate at 0.3 sec, corresponding to the max 
decceleration (1.75 gee’s), and it changes direction for the next 0.5 
sec until the value of 0.13 rad/sec is achieved. From this point the 
rate is decreased again until the constant value of -0.01 rad/sec 15 
achieved at 1.8 sec settling time . 


(2) The big overshoot of the rate is due to the produced pitch moment 
caused by the deflection of the damaged control surface while the 
other three are at their steady position. 


C. Pitch tail incidence (8, or DP), Figure 3.20, rad vs. time 
which has: 


(1) To achieve the acceleration response and due to the fact that 5, is 
locked at a negative angle of -0.2 rad, initally more control surface 
deflection is required (about 0.16 rad). At 0.3 sec, when the 
acceleration begins to decrease and the angular rate is zero, the tail 
incidence has the positive angle of 0.01 rad. A fluctuation begins 
before the deflection being stabilized at a constant angle of 0.06 
rad. 
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(2) Settling time equal to 1.8 sec. 
d. 08118000179 (61, 62, 635 64 or D's ), Figure 3.21, rad vs. time. 
This Figure illustrates the time response of the deflection 
of the four control surfaces, when one of them (5,) is locked at a negative 
angle of -0.2 rad, in order to achieve the acceleration response. Initially 
and when the 55, 63 and 54 control surfaces have zero angle of deflection 
ô, is already at the -0.2 angle. Then 55 and 5, are deflected at a negative 
-0.09 rad angle and 8; at a positive 0.09 rad. 
A comparison between Figure 3.21 and 3.1? indicates that when 
the surface ë, is hard over at a negative angle: 
(1) The undamaged control surfaces in order to counteract the failure 


and achieve the commanded acceleration, are deflected at a smaller 
angle. 


(2) The undamaged control surfaces are commanded to turn at angles 
with opposite sign comparing with those in Figure 3.17 ($, 
positively hard over). 


2. Damage One Control Surface 

Тһе ô; control surface was selected to simulate the inoperative 
actuator. The corresponding state variable 6, was therefore set equal to 
zero. The step input of one gee remained the same. The initial condition on 
the control surface was set at zero. This type of damage could also 
represent the distruction of the control surface. 

For the above case the time response plots of the state variables 
of interest were obtained and analyzed according to the defined 


performanse criteria in Section (I1I.B.4). Later, in Section 5, a comparison 
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Figure 3.14 Pitch Normal Acceleration vs. Time; Hard Over 6, Control 
Surface; 85508 5180 Іприї: п,.= 1988,5 1с=0.2 Гай, 62, ӧзо 5420. 
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Figure 3.15 Pitch Angular Rate vs.Time; Hard Over &, Control Surface; 
2.20; Step Input: n,,=1 gee, §,,=0.2 rad, Sn, 636, $4с=0. 
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Figure 3.16 Pitch Tail Incidence vs. Time; Hard Over 6, Control Surface; 
3,70; Step Input: 1051 068, 610-09.2 гад, 626 636 64:70. 
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Figure 3.1? Control Surfaces Deflections vs. Time; Hard Over &, Control 
Surface; 3,20; 5180 11001172051 098, 610-0.2 180, 620 630 64050. 
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Figure 3.18 Pitch Normal Acceleration vs. Time: Hard Over 64 Control 
Surface; a,=0; Step Input: n,.=1 gee, §,=-0.2 rad, 85, 53,8420. | 
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Figure 3.20 Pitch Tail Incidence vs. Time; Hard Over ô; Control Surface; 
35-0; Step Input: n, cz 1 gee, 64:2-0.2 rad, бос, 536, 20 
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Figure 3.21 Control Surfaces Deflections vs. Time; Hard Over ô; Control 
Surface; a,20; Step Input: n; cz 1 gee, 6,:2-0.2 rad, 62, 63, 6470. 
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between the different types of damages and the unimpaired system was 


performed. 


а. Pitch normal acceleration (nz or Nz), Figure 3.22, g’s vs. | 


time, which has: 

(1) Time constant or rise time 0.41 seconds. 

(2) A 0.005% steady state error. 

(3) Minor overshoot and delay time 0.35 sec and 

(4) A 1.6 sec settling time at which the acceleration response is 

achieved. 
b. Pitch angular rate (q or Q), Figure 3.23, rad/sec ys. time, 

which has: 


(1) A maximum body angular rate of 0.17 rad/sec is achieved at 0.3 
second. 


(2) Minor overshoot and a settling time equal to 1.6 after which the 
angular rate remains constant at -0.01 rad/sec. This is due to the 
fact that the transfer function 0/62 (Equation III.B.3-1) is a type 
zero system and has a finite error to a step input. 


с. Pitch tail incidence (6p or DP), Figure 3.24, rad vs. time, 
which has: 
(1) A 0.054 rad command from the actuator network. 
(2) А 0.4 вес delay time and negligible overshoot. 
0 08118000175 (84, 82, 85, 84 or D's), Figure 3.25, rad vs. time, 
which has: 
(1) A positive deflection angle of 0.035 rad for the & and 84 actuators 


while the 5: is deflected at a negative and equal angle 
(-0.035 rad), and $4 as damaged remains constantly at zero. 
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(2) A 1.6 sec settling time for each of the unimpaired surfaces. 


(3) A rise time equal to 0.7 sec and minor overshoot correspondingly 
for the undmaged surfaces. 


(4) The deflection angle of the unimpaired control surfaces has a 
pitching moment effect on the system. 
A comparison between the deflection angle in Figures 3.13, 3.17, 
3.21, and 3.25 indicates that for the same step input: 
(1) When one control surface is inoperative or shot away the system 


commands the undamaged surfaces at larger angles than those 
obtained with the unimpaired system. 


(2) When one control surface is inoperative or shot away the system 
commands the undamaged surfaces at smaller angles (in comparison 
with the hard over type of damage). 


3. Damage Two Control Surfaces 

In this case two control surfaces, on the same side of the 
missile were considered damaged (8, and 85). The two state variables 5, 
and 8, were therefore set equal to zero. The step input of one gee 
remained the same. The initial condition on the control surfaces was set 
at zero. Then the time response plots of the state variables were obtained 
and analyzed. 

8. Pitch normal acceleration, Figure 3.26, which has a 0.6 rise 
time, 0.38 sec delay time and 0.005% steady state error. Note the faulty 
behavior at the start of the response (0 to 0.15 sec). Because the system 
is a nonminimum phase system (i.e., a zero lies in the right-half s plane) 
the transient response starts out in the opposite direction to the input 


i.e., there is an inherent delay in the transient response of the system. 


69 








b. Pitch angular rate, Figure 3.27, which has a maximum rate 
of -0.16 rad/sec at 0.5 sec. Then increases until the constant negative 
rate (-0.01 rad/sec) is achieved at 1.3 seconds. 

c. Pitch tail Incidence, Figure 3.28, which has a 0.8 sec rise 
time and 0.3 sec delay time. The command from the actuator network is 
0.054 rad which is achieved at 1.6 sec settling time. 

d. 08116800175 (61, 55, 53, 54 or D’s), Figure 3.29, in which 
the two damaged control surfaces (5, and 55) remain always at zero 
(damaged) and the two others opposite deflected at the angle of 0.053 rad 
angle after a lapse time of 1.6 sec. Note that for this system with two 
impaired control surfaces more angle of deflection is required (in 
comparison with the system having one actuator shot away). 

4. Damage Three Control surfaces 

In the last case three of the control surfaces 4;, 45 and 53 were 
damaged or inoperative. 

With a step input of one gee’s for the commanded normal 
acceleration as in the unimpaired and other type damaged models, the 
plots of the state variables of interest were obtained and analyzed. 

Figures 3.30 through 3.33 illustrate the response of the output 
variables. Figure 3.50 shows that the commanded one gee acceleration can 
still be achieved in the absense of three control surfaces (84, 62, 63). 
However, some oscilation is evident into the response. The loss of the 
surfaces has produced a more lightly damped system as evidenced in 
Figures 3.50 through 3.55. The oscillatory and convergent behavior of the 
pitch rate and the undamaged surface (54), shown in Figures 3.51 and 5.55, 
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illustrate how the surface attempt to “fill in” for the absence of the 
inoperative control surfaces. specifically the response of the pitch normal 
acceleration (Figure 3.30) has: | 
(1) A 1.01 sec time constant. 
(2) A 0.38 sec delay time and 
(3) A 3.0 sec settling time at which the steady-state value of the 
commanded acceleration is achieved. 

A comparison with the response of the corresponding deflection 
of the unimpaired system (Figure 3.13) and the other inoperative cases, 
Figures 3.21, 3.25 and 3.29, indicates that as the number of the 
inoperative control surfaces increases, the deflection of the undamaged 
actuator(ís) also increases and more time is required to achieve the 
steady-state value. 

9. Comparison with the Unimpaired System 

In this section a performance comparison of each type of 
impairment with the undamaged system is performed and an analysis of 
their corresponding time response plots is presented. 

In order to use the N.P.5. EASYPLOT program the state variables 
had to be input as data. This was accomplished by writing a computer 
program ( Appendix D). 

with the above data Figure 5.54 is obtained. This corresponds to 
pitch normal acceleration (n, or NZ) for the unimpaired, hard over 8,, and 
inoperative or shoot away one or more control surfaces, configurations. 

The hard over case shows greater excursion in the normal 


acceleration with a peak value of -1.25 gee’s at 0.2 sec being achieved. 
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Figure 3.22 Pitch Normal Acceleration vs. Time; 
Inoperative &, Control Surface; 0] 
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Figure 3.23 Pitch Angular Rate vs Time; 
Inoperative &, Control Surface; a,=9; 


Step input: nzc =1 g8e, 610 2 0.2, 620 630 640 7 0 
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Figure 3.26 Pitch Normal Acceleration vs. Time; 
Inoperative 6; and 62 Control Surfaces; a,=0; 


step input: nze =1 gee, 84, = 0.2, 625, ӧз бас = 0 Габ. 
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Inoperative 61 and 52 Control Surfaces; a,=0; 
Step input: nz, =1 gee, бү, = 0.2, бәс, бз, бас = Or rad. 
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Figure 3.33 Control Surfaces Deflections vs. Time; 
Inoperative 8,, 8; and $3 Control Surfaces; a,=0; 
Step input: nz, =1 gee, $4, = -0.2, Son, 635, 640 = Ô rad. 
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The inoperative cases are all quite similar except for the last 
case where only one control surface is available. In this latter case an 
oscillatory behavior is evident. The performance criteria to the one gee 
command are summarized in Table 5.2. 


TABLE 3.2 


PITCH NORMAL ACCELERATION PERFORMANCE 


Unimpaired- Hard Over- Inoperative Control Surfaces 
One Two Three 


Time Constant(sec): 0.50 ....... 9 0.60....... 94600952০22 1.01 
Delay Time(sec): 0.28... Um 055.....0 5644 0.40 
Overshoot(3): 22% 0.01%.....0.01%....... 0.01% 


Figure 3.35 compares the unimpaired system with the three 
cases of control surfaces impairments. The curves include the following 
four cases: 

(1) Original system 

(2) Inoperative one control surface 
(3) Inoperative two control surfaces 
(4) Inoperative three control surfaces 

This Figure shows that as the number of the damaged control 
surfaces increases the delay time increases also from 0.28 sec for the 
unimpaired system to 0.41 sec for the system with three inoperative 


Control surfaces. The increase also of the time constant from 0.5 sec 
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(unimpaired system) to about 1 sec (three inoperative actuators), 
indicates that as the number of the damaged control surfaces increases 
more time is required to achieve the steady-state value of the normal 
acceleration. Also observe the initially faulty response of the nonminimum 
phase system to the step input which decreases as the number of the 
inoperative control surfaces increases. 

Figure 3.36, shows that to achieve the acceleration response the 
system with the hard over &, type of damage, requires more angular rate 
(q or Q), compared with the unimpaired and the other three damaged 
system cases. The produced pitching moment is evident. Initially and 
because the 6, control surface is deflected at the 0.2 rad angle while the 
other three are at the zero position, the first tenth of the second the rate 
increases from zero to 0.45 rad/sec. The rate then begins to decrease 
very rapidily to a rate of -0.38 rad/sec at 0.4 sec. 

Figure 5.5? compares the unimpaired system with three cases of 
control surfaces impairments. The curves include the following four cases: 

(1) Original system 

(2) Inoperative one control surface 
(3) Inoperative two control surfaces 
(4) Inoperative three control surfaces 

This Fígure shows that as the number of the damaged control 
surfaces increases, the angular rate of change is reduced, and more 
pronounced fluctuations occur, i.e., the system becomes more oscillatory. 
These fluctuations are connected with the pitching moment induced in the 


missile by the failed actuators. 


85 


Finally Figures 5.58 and 5.59 compare the unimpaired system 
with four cases of control surfaces impairments. The curves include the 
following five cases: 

(1) Original system 

)2( Hard over one control surface 

(3) Inoperative one control surface 
)4( Inoperative two control surfaces 
(5) Inoperative three control surfaces 

Figure 3.38 and 3.39 show that the required (steady-state) pifch 
tail incidence for the uncoupled longitudinal autopilot is 0.054 rad 


(unimpaired or damaged). 
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Figure 3.34 Comparison of the Pitch Normal Acceleration between the 
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Unimpaired System and the Four Cases of Impairment. 
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Figure 3.35 Pitch Normal Acceleration; Uncoupled open loop Channel: 
Unimpaired and Damaged Control Surfaces. 
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Figure 3.36 Comparison of the Pitch Angular Rate between the 
Unimpaired System and the Four Cases of Impairment. 
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Figure 5.5? Pitch Angular Rate; Uncoupled open loop Channel: 
Unimpaired and Damaged Control Surfaces. 
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Figure 3.38 Comparison of the Pitch Tail Incidence between the 
Unimpaired System and the Four Cases of Impairment. 
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Figure 3.39 Pitch Tail Incidence; Uncoupled open loop Channel: 


Unimpaired and Damaged Control Surfaces. 
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IV. LATERAL MODEL METHODOLOGY 


A. INTRODUCTION 
In this chapter the design and analysis of the continuous uncoupled 
lateral (i.e., yaw and roll) channels of a B.T T. autopilot derived in (Ref. 1 


and 4] is reviewed. À modified lateral autopilot is then designed with the 


. four control surfaces split into independent components. The model of 


elliptical (yaw channel) and circular (roll channels) configurations are both 
shown in Figure 4.1. Both airframes have the same cross sectional area 
distribution and specifically, the circular cross sectional body has a closure 
ratio Abase/Amax of 0.69 with Amax occuring at 68% missile body, whereas 
the elliptical airframe has a 3:1 cross section. 

A comparison between the two systems (original and modified) shows 
that they are identical. The transient response plots and their analysis 
demonstrate the performance of the original and modified yaw and roll 
channels. 

Having established the validity of the modified system one then uses 
the multiple-input/multiple-output control theory and state-space 
representation [Ref.3] to investigate the performance of the modified 
lateral autopilot for various classes of damages. 

For the investigation of the lateral channel, two classes of 
impairments at the control surfaces are considered: 

(1) Locked or hard over and 


(2) inoperative with the control surfaces centered or shot away. 
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The plots of the unimpaired modified yaw and roll autopilots are used 
as a basis for comparison. The impairment induced effects on roll angle, 


rate and tail incidence are assessed. 


B. UNCOUPLED LINEAR YAW CHANNEL (ORIGINAL) 
1. Aerodunamic Model, Airframe Configuration ang Control Lor 
A general block diagram of a B.T.T. autopilot with all its channels 
is shown in Figure 3.2. Inertial acceleration commands are applied in polar 
coordinates, i.e., magnitude of the command (nc) applied to pitch autopilot 
and the direction (9c) is applied to the roll autopilot. The yaw autopilot is 
slaved to the roll autopilot in order to minimize sideslip angle by 


coordinating the missile yaw and roll motion. Achieved maneuver plane 
accelerations in rectangular coordinates, (nz) and (ny) were determined by 


resolving achieved body-fixed accelerations, (т) апа (гу) through the 
missile roll angle (9), i.e., the Euler angles © and W are assumed to be 
sufficient small, so that the sideslip angle (6) is minimized and the motion 
between the yaw and roll channels is coordinated. 

A general diagram of the uncoupled yaw channel where the 
aerodynamic and yaw control law are involved is shown in Figure 4.2. In 
Figure 43 (Ref. 1] the yaw control law is presented. The yaw normal 
acceleration (yc) is not used to command the coupled B.T.T. autopilot, 
although, it is used for the design and analysis of the uncoupled channel. The 


command used by the coupled system is shown in dashed lines and is a yaw 
angular rate command (r,). The yaw control law is governed by missile body 
yaw angular rate (r) and yaw normal acceleration ( ny ). 
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The flight conditions are the same as was mentioned for the pitch 
channel, (i.e., 60 Kft altitude and Mack number 3.95). The yaw control law 
determines the required yaw tail incidence (&,,) command to the actuator 
network which is approximated as a first-order lag at 30 Hz. 

2. Requirements for the Classica! Design 


Requirements for the classical design approach are: 


a. High frequency attenuation in actuator command branch 2 15 dB at 
100 rad/sec and zerb angles of attack and sideslip. 


b. Relative stability: gain margins 2 6 dB, phase margins 230 deg with 
a goal of 12 dB and 50 degrees. 


c. Acceleration time response: 


(1) 63% time constant of approximately 0.4 sec for a step command of 
acceleration (ry) at flight conditions of interest and at a, and 8-0. 


(2) Overshoot « 108 for better sideslip control. 
(3) Steady state error not necessarily zero. 
3. Iransfer Functions of Aerodynamic Model and Analysis 
The aerodynamic transfer functions of the uncoupled yaw channel 
autopilot are according to Karadimas, C. thesis [Ref. 4]: 
a. Transfer function of yaw angular rate: 
| 











0.15 { $ + 1} 
1 0.0748 
— -—— (deg/sec/deg) /(ІУ.В.5-1) 
Sy I | 
( 6 + 4 $ + 1} 
4.3 -4.25 


b. Transfer function of Yaw Normal Acceleration: 
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l 1 











| 6% | () S+ i} 
Ty 11.48 -11.48 
z= —— O (q's/deg) (IV.B.3-2) 
Se 
{ S EM 
4.5 -4.25 


Rearranging the above transfer functions in terms of the variable 
pairs (r,8,) and (тё) respectivelu and applying the inverse Laplace 
transformation a set of linear differential equations can be obtained. 
Then, using rules of state-space representation, the following equations 


are obtained [Ref. 4]: 


X, = X, - 29.2321 6, (1У.В.3-3) 
כא‎ = -16.7224% | - 0.0434X, + 0.26875, (IV.B.3-4) 
Z1 = 2; - 0.08018, (IV.B.3-5) 
Zo = -16.7224, - 0.04352, - 150.77248, (IV.B.3-6) 
Ү = 3.35742 -4Y * 6.43418, - 3.5574 (IV.B.3-7) 
$Yc- 6.06X, + 0.606X, + 2.04132, + 3.648Y - 13.8612 

-2.0413те (IV.B.3-8) 
SY = 188.48,, - 188.46, (IV.B.3-9) 


Using state-space representation, the Equations (IV.B.3-3) through 


(IV.B.3-9) can be modeled in a seventh-order system of the form: 


X = FX + GU 
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with the continuous plant system F and control matrix G as follows: 


X, 0. j 0 0 0 0 -29.2321] ( X, 
X» -16.722 -0.0434 0 0 0 0 0.26870] | x, 
21 0 0 0 | 0 0 -0.08010| | 2, 
| = 0 0  -16.?22 -0.0435 0 0 -150.772 | | 2, 
Y 0 0 3397 «0 -4 о 6.434 || Y 
86| | 6.06 0.606 2.041 0 3.646 0  -13.86 | | 
[ww—Fli“ оо 0 0 о 188.4 -188.4] | 6] 
| 9 1 (те 
|| 
0 
. + 0 
- 3.3574 
- 2.0413 
0 


where the state variables are: 
- у, Хә : yaw angular rates 
- Z;, 2 : yaw normal acceleration 


ary : output of acceleration compansator 
- 8৬০ : input command in the actuator 


- бу : yaw tail incidence 


Using a step input function which represents the `] gee command 
at ае = 0, the time response plots (Figure 44), the Yaw Normal 
Accelerarion (ny) (Figure 4.5), and the Yaw Angular Rate (В) and the Yaw 
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Tail Incidence (Sy) (in Figure 4.6) were obtained. The transient response 


satisfies the performance criteria and definitions given at Section (III.B.4), 
and are identical with the results presented in [Ref. 4]. 
In Figure 4.4 the Yaw Normal Acceleration is plotted which has: 
(1) A 0.31 sec time constant. 
(2) À 0.3 sec delay time and 6.4€ overshoot. 
(3) À 0.018 steady-state error. 
(4) Note the faulty behavior at the start of the response (0 {0 0.2 sec). 
Because the system is a nominimum phase sustem (i.e., a zero lies 
in the right-half s plane) the transient response starts out in the 


oppsite direction to the input i.e., this is a inherent delay in the 
transient response of the system. 


In Figure 4.5 the Yaw Angular Rate is presented which has: 
(1) A maximum body angular rate Of 0.6 rad/sec achieved at 0.3 second. 
(2) А minor undershoot and 2.2 sec settling time after which the 
angular rate remains constant at 0.01 rad/sec. This is due to the 


fact that the transfer function r/&y (Equation IV.B.3-1) is a type 
zero system and has a finite error to a step input. 


Finally Figure 4.6 shows that a negative angle of deflection of 
-0.14 rad is required in order to achieve the commanded step input 
acceleration. 

The above analysis has revealed that the performance of the system 
confirms to the previous mentioned requirements. 

From the pole-zero plot in Figure 2.? [Ref.4] we can see that the 
continuous open loop system is stable, since the s-plane poles are: 

Sı = -175.010 -jO (IV.B.3- 10) 


98 


99 





(91095 9/1) 501124161 ШшШ027 59012011115 12213101113 pue 101104117 10 9190011 l't 3161 


£69 ари P 
: 1797 Арос шпинхеиј וט וטכןו59‎ 
| OCZ SUOISUALUIP үу 





/ bee 


„ОЕ 





100 





Do 
“је. 7% 
+ 
ДЕ 
: E 
(бар) .. | £ 4G 
В 5 
+ 
4 
А 
A כ‎ sb 
(бэр)^ $ = 
Є`{6 | 


(ре!) ме| )5.5| א‎ 
олио? 
"AQ MP A 





4910110 


(2а5/рез) | 
| 


101 


Alt 


А01 10213007 АВА <р Ә1ПВІ4 


12511913 
48 [D2112) 


QUIC) JI V 


(2as/pe1)1 








0 


9рбІСО 





0. 


Autopilot; Continuous Open Loop system; Elliptical Airframe; 
102 


Figure 4.4 Yaw Normal Acceleration vs. Time; Uncoupled Yaw Channel 
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Figure 4.5 Yaw Angular Rate vs. Time; Uncoupled Yaw Channel 
Autopilot; Continuous Open Loop System; Elliptical Airtrame; 


Step Input HY = 1 gee; a, = 9. 
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Figure 4.6 Yaw Tail Incidence vs. Time; Uncoupled Yaw Channel 
Autopilot; Continuous Open Loop system; Elliptical Airframe; 


Step Input HY = 1 gee; a, = 9. 
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5) = -5.7665 4 (IV.B.3-11) 


5ҙ < -5.7665 -|10.4804 (IV.B.3-12) 
S4 = -0.0217473 *j4.08919 (IV.B.3-13) 
Ss = -0.0217473 -j4.08919 (IV.B.3-14) 
56 = -2.9484 +]2.99489 (IV.B.3-15) 


57 = -2.9484 -]2.99489 )1/.8.3-16( 


C. UNCOUPLED LINEAR OPEN LOOP ROLL CHANNEL (ORIGINAL) 


The Aerodynamic Model, the Airframe Configuration and the‏ .ו 


Control Law 

Generally, the purpose of a roll channel autopilot is to command 
a roll ($) signal, in order to roll the missile through an angle (9) 
measured from the vertical axís. | 

While the pitch channel acceleration is commanded to produce 
the total magnitude of the guidance acceleration command, the roll 
channel of the B.T.T. autopilot is commanded to roll the missile so as to 
put the prefered maneuver direction in the direction of the guidance 
acceleration command, as is shown in Figure 3.2. 

À block diagram of the uncoupled roll channel is shown in Figure 
4.7, and the roll control law in Figure 4.8 ( both Figures taken from 
[Ref. 1]. In these diagrams the same aerodynamic model is used. 

As is shown in Figures 4.7? and 4.8 the control law is commanded 
by a roll angle (c), is governed by roll angular rate (P) and roll angle (9) 
and determines the required command (&,) to an actuator. This is 


approximated as a first order lag at 188.4 rad/sec. 
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B 


The circular airframe configuration is chosen for this channel 


and is shown in Figure 5.1. The flight conditions are the same as for the . 


Yaw and Pitch channels, i.e. 60000 ft and Mack number 3.95. 
2. Requirements for the Classical Design 
Requirements for the classical design approach are: 
а. High frequency attenuation in Actuator command branch » 15 dB at 


100 rad/sec and zero angle-of-attack. This requirement limits the 
speed of roll angle response. 


b. Relative stability: gain margins 2 30° with a goal of 12 dB and 50°, 
i.e. the same as for the yaw channel. 
c. Acceleration time response : 


(1) 63% time constant 0) 0.5 seconds for a step command of roll 
angle at the flight condition of interest and zero angle of attack. 


(2) Overshoot < 10%. This requirement is conservative and may be 
relaxed depending upon guidance results. - 


(3) Zero steady state roll angle error. If not zero, it will affect 
guidance. 


3. ransfer Funcions of the Aerodunamic Model and Ane 
The aerodynamic transfer function of the uncoupled roll channel 
autopilot are (Ref. 4]. 
a. Transfer functions of the Roll Angular rate: 
P=(qSd/1L,)*57.3*C)m Š (rad/sec?) 


where: 

4 - 1650 (Ibs/ft?) 
9 - 3.14159 ft? 
Ч = 21% 
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> =40 
Cy m =0.031 
Р = 8.0346 бр (IV.C.3- 1) 


b. Transfer function of Roll Angle: 
ÓP = Р (rad/sec) )1/.6.3-2( 
Rearranging the above equations along with those given in 


(Ref. 1] for the actuator and using rules of state-space representation, the 


following equations are obtained [Ref. 4] and verified: 


P = 8.034628, | (I¥.C.3-3) 
j =P | (IM.C.3-4) 
Х = -17.6ф -8Х +17.6%, (IV.C.3-5) 
Y = -0.75325P -0.88359 +0.35165X -5Y -0.403346, 

+0.883529, (IV.C.3-6) 
X, = -6X, +0.62622918, (IV.C.3-7) 
bp, = -0.1027 1P -0.120489 +0.04795X +14.3182Y 

-14.18184*[ - )1/.6.3-8( 
ӛн - 10795,326р, - 188.46; (IV.C.3-9) 


Utilizing state-state representation, the above Equations can be 
modeled as a seventh-order system in the format X - FX *6U with the 
continuous plant system F and control matrix G as shown in Table 4.1 and 


the state vector with the corresponding state variables as follows. 
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Р Roll angular rate 

9 Roll angle 

X Output of roll angle compansator 
Y Output of roll rate compansator 
Ay Ouput of pseudodifferensiator 
бре Input command in the actuator 


бр Roll tail incidence 


Executing the interactive implementation of the optimal systems 
control design program (OPTSYSX) using an input step function which 
represents “1 radian command” at zero trim angle of attack, the graphical 
time response plots were obtained for the state variables of interest. 

Figure 4.9 shows the Roll Angular Rate (P) (rad/sec vs. time), 
Figure 4.10 the Roll Angle (9) (rad vs. time), and Figure 4.11 the Roll Tail 
Incidence (85) (deg. vs. time). A close observation in Figure 4.10 indicates 
that for the commanded input vector the time response plot of the roll 
angle has a 0.55 sec time constant, 5.3% overshoot and a steady-state 
error equal to zero. These results are in accordance with the requirements 
mentioned in the previous section. 

From the pole-zero plot of Figure 2.23 [Ref.4] we can mention 
here that the continuous open loop system is stable, since the s-plane 


poles are: 


S4 = -174.785 -j0 (IV.C.3-10) 
<) = -9.2510 8 (IV.C.3-11) 
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$3 = -9.25097 -]28.4098 (IV.C.3-12) 


54 = -2.46608 +2 (IV.C.3-13) 
$5 = -2.46608 +j2.71152 (IV.C.3-14) 
56 = -8.98209 +10 (1У.Г.5-15) 
$7 = -5.19932 +10 (IV.C.3-16) 


4. Design Approach and Analysis of the Modified Model (Unimpaired) 


To study reconfigurable flight controls the missile model should 
be capable of using the set of four available control surfaces to produce the 
required forces and moments. | 

То achieve this each control surface (84, 85, 65, 64) was considered 
as an independent component. | 

The equations were then coupled so that a component of control 
surfaces used alone or in conjunction with another control surface, will 
produce the appropriate moments. 

Taking under consideration the conventional control surface theory 
of Chapter II. for the Sank-to-Turn missiles, and for the circular 
configuration shown in Figure 3.4, one can easily show, for the continuous 
open loop roll autopilot that the input command in the actuator network 
(Spe) can be split up into independent components, using the following 
equation which computes the deflection inputs to Equation (IV.C.5-8). Also 
refer to Figures 2.2 and 3.4 which show the control surface position and 


corresponding positive deflections. 


+ - CT = Š2 = 63 *64) / 4 (IV.C.4-1) 
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Replacing the input vector in Equation (IV.C.3-8) and (IV.C.3-9) 


with the above Equation and applying the superposition principle, gives the 
following modified airframe equations: 


бр = -160.45р + 10795.32 (5. - 6 - 63 + 64) /4 (IV.C.4-2) 


(61 -62 -63 +84) / 4 2 -0.1027 1P -0.120489 +0.04795X +14.318Y 
-141818X, -15/4 (84+ב82-83-,8)‎ - 7 
- 0.12048%- )1/.6.4-3( 


In Figure 4.6 the modification of the commanded tail incidence is 
presented. 
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Figure 4.6 (a) Original System; (b) Modified System. 


Rearranging the above equations and using the transfer functions 


of the system for the original uncoupled roll channel, Equations (IV.B.3-3) 
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through (IV.B.3-9), the following Equations were derived and can be 


modeled in a ten order system for the rol! channel: 


(IV.C4-4)‏ = ןק 
j =P )1/.6.4-5(‏ 
Х = -17.6ф -8Х +17.64г )1/.6.4-6(‏ 
Ү --0.75325Р -0.88352% %0.35165Х -5Ү‏ 

-0.403345p «40.883529 c )1/.6.4-7( 
১1 = -6%| 8% (1V.C.4-8) 
$, = -0.10271Р -0.120489 +0.04795Х + 14.5182% 

-14.18184Х, -158, -0.140628р +0.120489е (IV.C.4-9) 


$2 = +0.10271Р +0.12048$ -0.04795Х - 14.3182 
+14.18184Х; +155, +0.140625р -0.12048. (IV.C.4-10) 


$, = +0.10271Р %0.12048% -0.04795Х - 14.3182У 
+14.18184Х +155, +0.140625р -0.12048%г  (IV.C.4-11) 


4 = -0.1027 1P -0.120489 +0.04795X + 14.3182Y 
-14.18184X, -156, -0.140625p %0.12048%- (IV.C.4- 12) 


5ң - 2698.835) -2698.835; -2698.835; %2698.835, 
-188.46ң (IV.C.4-12) 


We thus have a multiple-input-multiple-output (MIMO) system 
with the following order: 
- X: 10x1, the time derivative of X matrix. 
-F:10x10, the open loop dynamic matrix or plant. 
- X: 10x1, the variables assignments matrix. 
- G: 10x5, the control distribution matrix. 


- ש‎ : 5x1, the control input matrix. 


112 


The above control and output matrices of the Uncoupled Modified 
Unimpaired Roll Channel Autopilot are presented at the Appendix B. 

The Eigenvalues found executing the OPTSYS Control Program, 
that is the DET(sI-F) values, show that the open loop system is stable, 


since the s-plane poles are: 


91 = -185.057 + jo (roll angular rate) 
5, = -9.27462 + j29.2679 (roli angle) 
בָ5‎ = -9.27462 - 9 
54 = -2.3582! + 6 
95 = -2320082] - 6 
56 = -8.89209 + ]0 (input command in the actuator) 
57 = 18517.ם-‎ + 0 | (roll tail incidence) 
Зв = -15.000 + jo 
53 = -15.000 + jo 
510= -15.000 + jo 
о. Analysis and Comparison with the Original System 
Applying the same step input function as in Section IV.C.3 which 
represents the 71 radian command’, and for zero rad deflection command 
for the four control surfaces (1.е, бух 625 63১ 64. - 0) the time response 
plots of the following state-variables of interest were obtained and 
presented in the Figure 4.12 through 4.15. 
a. Roll Angle (9 or F), Figure 4.12, rad vs. time which has: 
(1) A time constant 0.6 seconds 
(2) A 5.33 overshoot. 
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(3) Zero steady-state error and so the guidance will not be effected. 
(4) Delay time 0.48 sec and settling time 2.2 sec when the steady 
state roll angle is achieved. 
b. Roll Angular Rate (P or p), Figure 4.13, rad/sec vs. time, 
which has: 
(1) A maximum body angular rate 1.58 rad/sec at 0.4 second. 
(2) A 4.2% undershoot and a 2.4 sec settling time. 
c. Roll Tail Incidence (Sp or DR), Figure 4.14, degrees vs. time, 
which has: | 
(1) Initially and very fast, 0.1 sec, the maximum required deflection 


occurs, that is 1.1? degrees. Then changes direction to the opposite 
deflection angle in 0.6 sec to -0.52 deg. before stabilizing. 


(2) A 2.6 sec settling time . 
d. 06118000175 (5), 85, 85, 84 or D's), Figure 4.15, rad vs. time, 
which has: 
(1) The actuators 8, and 84 have the same sign (positive), as was 
expected according to the sign convention of Figure 2.2, and are 
first deflected to a 0.02 rad angle and then to -0.05 radians after 


0.6 seconds. Finally the steady state zero angle of deflection is 
achieved at 2.6 seconds. 


(2) The & and 53, have the same response sa above but with a negative 
sign. 
As was shown the relative stability for zero-angle-of-attack is 
satisfactory for the modified unimpaired system. 
A close observation of Figures 4.12 throught 4.15 and the above 
analysis, indicates that the corresponding obtained responses are 


identical, with those of the original roll channel, Figures 4.9, 4.10 and 
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4.11. So, the performance of the designed system meets the perviously 
mentioned requirements. 
In table 4.2 a summary of the stability margins is presented, 


obtained from the frequency response plots calculated up to 100 rad/sec. 


TABLE 4.2 


PHASE CROSSOVER FREQUENCIES AND GAIN MARGINS; 
UNCOUPLED MODIFIED ROLL CHANNEL(UNIMPAIRED); 
CIRCULAR AIRFRAME 


PHASE CROSSOVER GAIN 
EREQUENCY(RAD/SEC) MARGIN(DB) 
ROLL ANGLE (9) 10.0461 60.2165 
ROLL ANGULAR RATE (P) 10.0461 60.2165 
| ROLL TAIL INCIDENCE (8р) 10.0461 42.1172 


We procceed now to the analysis of the impaired control surfaces 
for the modified system. 


D. CONTROL EFFECTIVENESS DUE TO THE CONTROL SURFACES 
IMPAIRMENT AND COMPARISON 


The time response plots for the system are obtained in this section 
for two types of damages. The first type is the locked or hard over 
actuator at a small angle and the second type includes the centered or 
inoperative one, two and three control surfaces. 

Note that the first phase of the design approash for the CBTT 


autopilot was to design each channel independently with all coupling 
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Figure 4.9 Roll Angular Rate vs. Time; Uncoupled Roll Channel; 
Continuous Open Loop System; Circular Airframe; 


Step Input $, = 1 rad; a, = 0. 
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Figure 4.10 Roll Angle vs. Time; Uncoupled Roll Channel; 
Continuous Open Loop System; Circular Airframe; 
Step Input f, = 1 rad; a, = Q. 
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Figure 4.11 Roll Tail Incidence vs. Time; Uncoupled Roll Channel; 
Continuous Open Loop System; Circular Airframe; 


Step Input $0 - ! rad; a, = 0. 
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Figure 4.12 Roll Angle vs. Time; Uncoupled Modified Unimpaired Channel; 
Continuous Open Loop system; Step Input : 
$e = 1 and 616 62e 636 646 = 0 , дв = 0. 
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Figure 4.15 Roll Angular Rate vs. Time; Uncoupled Modified Unimpaired 
Channel; Continuous Open Loop system; Step Input : 
Ус = | and 910 620 630 бас = 0 ? 86 = 0. 
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Channel; Continuous Open Loop system; Step Input : 
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Figure 4.15 Control Surfaces Deflections vs. Time; Uncoupled Modified 
Unimpaired Channel; Continuous Open Loop system; 
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between channels removed. Sufficient high frequency attennuation was 
added for actuator and missile design so that the resulting missile body 
angular rates and control surfaces motion would represent a practical 
missile design. A relationship was established among the relative speeds 
of response of the uncoupled channels in order to meet CBTT requirements. 
The cross-coupling due to control surfaces impairment between the 
longitudinal and lateral motion of the missile is ignored when each control 
channel is treated independently. Cross-coupling of the controls is thus 
considered in Chapter V. 
1. Hard Over One Control Surface 
One control surface (8,), was selected to be deflected or hard 
over, i.e. it remained at a fixed angle while the other three control 
surfaces had a variable deflection angle. 
Using a step input 1 radian for the commanded roll angle (97) (as 
in the original and unimpaired models); an angle of 0.2 rad for the 
damaged surface (5;,); and a zero angle of deflection, for the remaining 


other three (620 53,, andS4.). The time response plots of the state 
variables of interest were obtained and analyzed with the performance 
criteria as defined in section (III.B.4). 
a. Roll Angle ( or F), Figure 4.16, (radians vs. time) has: 
(1) A time constant 0.6 sec and 0.3 sec delay time . 
(2) A constant 0.4 rad steady state error. 
(3) A 5.8% overshoot at 1.2 sec. 


(4) The roll angle is stabilized at an incorrect angle of 1.4 rad instead 
of the commanded 1 radian due to the impairment of 6, surface. 
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b. Roll Angular Rate ( P ), Figure 4.1? (radians/sec vs. time) 


which has: 


(1) Initially and due to the fact that 5, is locked at a positive angle of 
0.2 rad, to achieve the the roll command, the angular rate increases 
very fast the first tenth of the second to a 2.1 rad/sec rate. Then, a 
pronounced oscillation occurs until the 0.4 sec when the rate begins 
to decrease crossing the zero rate at 1.2 sec corresponding to the 
max roll angle (1.5 rad). 


(2) The initial increase of the rate is due to the induced rolling moment 
caused by the deflection of the damaged control surface while the 
other three are at their steady position. 


(3) A 6.6% undershoot at 1.5 sec and settling time 2.5 sec. 
c. Roll Tail Incidence (6a or DR), Figure 4.18 (deg. vs. time) 


which has: 


(1) At time equal to 0.1 sec the actuator is commanded to an angle 4.75 
degrees. This is is due to the fact that 8, is hard over at the 0.2 
rad angle while the other three remain have a zero initial condition. 
Then the tail incidence is stabilized after a 2 second period during 
which a fluctuation occurs. 


(2) A zero steady-state error and 2.2 sec settling time. 
d. Control Surface Deflections (8,, 55, 53,54 or D’s), Figure 
4.19, rad vs. time, which has: 


(1) The actuator 8,, is deflected at the angle of 0.2 rad, from zero sec 
up to 3 seconds. 


(2) The commanded input vector now causes the two surfaces ô, and Š; 
to turn in a positive angle of 0.15 rad, instead of a negative angle 
as when no impairment existed. | 


(3) A 1.6 seconds settling time to achieve this angle. 
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(4) The fourth surface 54 is commanded to turn at an angle of equal 
amount but with an opposite sign, for the same reasons. 


A comparison between Figure 4.19 and 4.15 of the unimpaired 
system, shows that to achieve the roll response the impaired sustem 
requires more control surfaces deflection and more time. 

2. Damage One Control Surface 

The ô, control surface was selected to simulate the inoperative 
actuator. The corresponding state variable 8, was therefore set equal to 
zero. The step input of one radian remained the same. The initial condition 
on the control surface was set at zero. This type of damage could also 
represent the dedtruction of the control surface. 

For the above case the time response plots of the state variables 
of interest were obtained and analyzed. Later in Section 5, a comparison 
between the different type of damages and the unimpaired system is 
performed. | 

a. Roll Angle ( 9 or F ), Figure 4.20, rad vs time, which has: 

(1) A 0.6 second time constant. 

(2) An overshoot equal to 5.6% and 

(3) Zero Steady state error. 

(4) Delay time 0.5 sec and 2.1 sec settling time at which the the roll 
response is achieved. 

b. Roll Angular Rate ( P ), Figure 4.21, radian/sec vs. time, 
which has: 


)1( A maximum body angular rate 1.58 rad/sec at 0.4 sec. 
(2) A 6.53 undershoot. 
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(3) A 2.5 sec settling time at which the zero steady state value is 
achived. 


c. Roll Tail Incidence ( 5 or DR ), Figure 4.22, deg. vs. time, 
which has: 
(1) At the very beginning of the response, 0.1 sec, the maximum 


deflection occurs, that is 1.1 degrees, to change then to the 
opposite deflection angle in 0.7 sec to -0.32 deg. 


(2) A zero steady state error and settling time 2.6 sec. 
d. Control surfaces deflection (Sm or DR), Figure 4.25, deg. vs. 
time which has: 


(1) The surface 5, as damaged remains constantly at zero angle of 
deflection. 


(2) The actuators 85 and ás are deflected at a negative angle 0.026 rad 
and are stabilized after 2.4 seconds. 


(3) The actuator 54 is deflected at an. equal but opposite angle as it was 
expected according to the selected sign convention. 
А comparison between the deflection angle in Figures 4.15, 4.19 
and 4.25 indicates that for the same step input: 
(1) When one contro! surface is inoperative or shot away the system 


commands the undamaged surfaces at larger angles (in comparison 
with the unimpaired system). 


(2) When one control surface is inoperative or shot away the system 
commands the undamaged surfaces at smaller angles than obtained 
with the hard over type of damage. 


(3) The required time that the undamaged control surfaces need to 


achieve the zero steady state value is the same (in comparison with 
the unimpaired system). 
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5 Tw ntrol Surface 


In this case two control surfaces, on the same side of the 
missile were considered damaged (8, and ô»). The two state variables ô; 
and 52 were therefore set equal to zero. The step input of one radian 
remained the same. The initial condition of the control surfaces was set 
at zero. Then the time response plots of the state variables were obtained 
and analyzed. | 
a. Roll Angle ($ or F), Figure 4.24, rad vs. time, which has: 
(1) Time constant 0.59 seconds 
(2) An overshoot 6.8% at 1.25 sec and an undershoot of 0.003% at 2.2 
seconds. 
(3) A zero steady-state error and a 0.5 sec delay time. 
(4) The commanded roll angle is achieved after 2.8 seconds. 
b. Roll Angular Rate (P), Figure 4.25, radian/sec vs. time, 


which has: 


(1) The peak body angular rate 1.61 rad/sec is achieved at 0.41 sec. 
(2) An 8.2% undershoot at 1.6 sec and 
(3) A 2.6 sec settling time. 
c. Roll Tail Incidence ( ác or DR ), Figure 4.26, deg. vs. time, 
which has: 
(1) The peak angle is 1.01 degrees at 0.1 sec and then in 0.6 seconds a 
negative angle (-0.35 deg) is commanded. 
(2) A second 6.3% overshoot at 1.6 sec and 


(3) Zero steady state error (settling time 2.6 sec). 
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d. Control surface deflection (84, 65, 63, 54 or.D's), Figure 4.27, 
rad vs. tíme, which has: 


(1) The two inoperative or shot away control surfaces 61 and 62 remain 
always at zero angle of deflection. 


(2) The two remaining actuators 6; and 54 moving at opposite 
directions, as expected, with a peak value of +0.35 radians. 


(3) The impairment of ô; and ô, has produced a more overshoot lags 
system as evidenced in Figure 4.27. 

(4) The convergent and oscillatory behavior of the two unimpaired 
control surfaces which attempt to "fill in^ for the damaged surfaces 
is shown in the same Figure. Also note that for this system with 
two impaired control surfaces more angle of deflection is required 
(incomparison with the system having one actuator shot away). 


4. Damage Three Control surfaces 

Іп the last case three of the control surfaces ô, ô> and ô; were 
damaged or inoperative. | 

With a step input of one radian for the commanded roll angle (92) 
as in the unimpaired and other type damaged models, the plots of the state 
variables of interest were obtained and analyzed. 

Figures 4.28 through 4.31 illustrate the response of the output 
variables. Figure 3.50 shows that the commanded one radian roll angle сап 
still be achieved in the absence of three control surfaces (6,, 62, 63). 
However, some oscilation is evident in the response. The loss of the 
surfaces has produced a more lightly damped system as evidenced in 
Figures 4.28 through 4.31. The oscillatory and convergent behavior of the 
roll rate and the undamaged surface (54), shown in Figures 4.29 and 4.31, 


illustrate how the surface attempt to ‘fill in’ for the absence of the 
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inoperative control surfaces. Specifically the response of the roll angle 
(Figure 4.28) has: 

(1) À 0.55 sec time constant. 

(2) A max. 10.5% overshoot at 1.2 sec and a 0.5% undersoot at 2.2 sec. 


(3) The commanded roll angle is achieved after 2.8 sec. 


(4) 80.51 sec delay time . 

À comparison with the response of the corresponding deflection of 
the unimpaired system (Figure 4.15) and the other inoperative cases, Figures 
4.19, 4.25 and 4.31, indicates that as the number of the inoperative control 
surfaces increases, the deflection of the undamaged actuator(s) also 
increases and more time is required to acheive th steady-state value. 

9. Comparison with the Unimpaired System 

In this section a performance comparison of each type of 
impairment with the unimpaired system, was performed and an analysis of 
their corresponding time plots is presented. 

For the above purposes, a computer program was written 
(Appendix D) to create a data file, for the state variables of interest, to 
be used as input data in the use of the EASYPLOT program of the NPS 
computer. 

For the decoupled rol] channel the first plots are for the roll 
angle and are presented in Figures 4.32 and 4.33 

In Figure 4.52 which compares the two types of impairments 
with the original (unimpaired) system it is quite obvious that the worst 


case is the hard over case because when the autopilot is commanded at a 
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roll angle ($c) equal to ! radian, the missile turns at a roll angle (9) of 
1.4 rad instead of 1 rad as for the unimpaired autopilot. 

Other characteristic point is that for the hard over case, the roll 
angle reaches its final value of 1.4 rad (80.29) in approximately 2 seconds | 


with a 5.8% overshoot, without any steady state error. 


| 
Figure 4.33 compares the roll angle of the unimpaired sustem 
with the three cases of control surfaces impairments. The curves include 
the following four cases: 
(1) Original sustem 
(2) Inoperative one control surface | 
(3) Inoperative two control surfaces 
(4) Inoperative three control surfaces 
This Figure shows that as the number of the damaged control 
. Surfaces increases the overshoot and the settling time icrease. | 


Figures 4.34 and 4.35 compare the roll angular rate of the 


unimpaired system with four cases of control surfaces impairments. The 
curves include the following five cases: x 
(1) Original sustem | 
(2) Hard over one control surface | 
(3) Inoperative one control surface 
(4) Inoperative two control surfaces 
(5) Inoperative three control surfaces 
In Figure 4.34 it is shown that the hard over case induce greater 
excursion in the roll rate with a peak value of 2.2 rad/sec at 0.38 sec 


being achieved. Due to the fact that 8, control surface is deflected at the 
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0.2 rad the initial rolling moment in the first tenth of the second 
increases from zero to 2.1 rad/sec. The other control surfaces have zero 
as an initial condition. 

Figure 4.55 compares the unimpaired system with three cases of 
control surfaces impairments (inoperative one, two and three actuators). 
In this is shown that the inoperative cases are all quite similar except for 
the last case where only one control surface is available. In this last case 
a larger delay time, overshoot and undershoot are evident. 

Finally in Figures 4.56 and 4.5? the roll tail incidence is 
compared and include the following cases: 

(1) Original system 

(2) Hard over one control surface 

(3) Inoperative one control surface 
(4) Inoperative two control surfaces 
(5) Inoperative three control surfaces 

Figure 4.56 shows that the required roll incidence for the hard 
over case is larger in comparison with the unimpaired system and with 
the other three cases of control surfaces impairment. In Figure 4.3? the 
maximum peak overshoot decreases with the number of control surfaces 


that are impaired. 


155 


1.6) 


.5 


э ө е ө өе е ч е е ө е є э аз ә е ® = = =,е е е = оо оо ооо ооо оо ооо о ооо ооо ооо ое е до а а = е = р = @ т е ө о ч ао оо оа оо = е э ® а оа е еэ е э оное. 


1 


1.4 


е ө е е е ® э = = е еэ = от = е = = е ® = ® = э эъ э а е сс о ооо ооо ооо оо ® = е ® = = = == = 


1.3 


-езетвейееесвес Жееееееейеевее-с фо ооо ооо ооо о Ѓѓе е оо одо е ев ово оо ое ао Бе э в еэ е дф э ө оо ое е ево со ово се е 


Pe 


6000800000 505 הש‎ ও ও ও ৪৮ со оо ое ও b ое 9» o 9 9» € 9 & e 9 9 9? 9? Ооо о ое а оо а ое то ом э ® ое о ое оо "9 da E от = = а e ое одо e o ^ 9? e e ма а е es 


1.1 


1.0 


0.9 


1 1 1 ' 1 1 1 ' 1 
. . . . , , ' . . П . . 1 . 
ооо обо ооо вое фео е 1.9.9 4... ....... ₪0. еее еве "еве ве ге k e... ... "donc anesson o еее о l wees... 
' ' . . . 1 1 1 ' . . 1 . . 
. . ' * ' . ' ' . . * . 0 
1 “ 1 . 1 1 4 


ооо ооо оо ще = = = ®1= ® = э э = = ® = 


0.8 


ооо оц ооо ооо ооо ооо ооо су э = = = эъ үе ооо о цело оо е ој оо зо ош ор оо ое и е е ци а а ое оо е о ооо оо оо ое е ор ое оо оо... 


F,ROLL ANGLE- RAD 
0.7 


.9..: 0 0. е фе чес је сос соло је е ое ло оо оа оо е оо оо 


= LEGEND 
ZJ- S  F,OROLL ANGLE-RAD 





0.0 0.4 0.8 1.2 1.6 2.0 2.4 2.8 
TIME-SEC 


Figure 4.16 Roll Angle vs. Time; Uncoupled Modified Roll Channel 
with 6, Hard Over; Continuous Open Loop System; 
Step Іприќ: ф. = 1 апа 5.20.2 rad; 62630645098 86 0. 
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Figure 4.1? Roll Angular Rate vs. Time; Uncoupled Modified Roll Channel 
with 5; Hard Over ; Continuous Open Loop System; 
Step Inputs: 6,;=0.2 rad; 64,53,54=0; Pe=1 rad and a, = 0. 
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Figure 4.18 Roll Tail Incidence vs. Time; Uncoupled Modified Roll Channel 
with 6, Hard Over ; Continuous Open Loop System; 
Step Inputs: $,20.2 rad; $2,852,420; $c-1 rad and a, = 0. 
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Figure 4.19 Control Surfaces deflections vs.time; Uncoupled Modified Roll 
Channel with 6, Hard Over ; Continuous Open Loop System; 
Step Inputs: 6,=0.2 rad; 65,53,54=0; Pe=1 rad and a, = 0. 


157 





0.6 0.7 0.8 
. 
. 
4 
. 
4 
0 
. 
. 
. 
, 
. 
. 
L] 
3 
. 
0 
. 
' 
4 
. 
ו‎ 
' 
. 
, 
. 
. 
* 
s 
, 
. 
4 
₪ 
4 
4 
r 
4 
. 
. 
. 
. 
0 
- 
. 
. 
4 
. 
. 
. 
t 
+ 
0 
. 
. 
. 
, 
0 
* 
0 
. 
. 
. 
, 
г 
0 
* 
* 
. 
* 
. 
P 
. 
. 
. 
. 
. 
. 
- 
. 
. 
. 
. 
. 
. 
ו‎ 
. 
. 
. 
, 
0 
+ 
е 
. 
. 
. 
. 
. 
. 
г 
, 
+ 
, 
* 
. 
. 


F,ROLL ANGLE- RAD 


0.5 


. . ' 
э ә ә 9. 9 9 9 а. 9 е э ое = э ое = ф =» е е ә еэ е ое э = ә =» э э = ор ое е ө == Ф = Фе ое оо ро е ое еэ а ое е ко е ооо оо ооо о ох ооо оо ооо оо ооо оо ооо ооо ооо оо оо ооо ооо 
1 1 1 


--е-еедйдеетее-е ”Ы--ә-с-есеібте «е-е «-ы" ооо ооо о ооо фо е-ес е одо оо ом оо оо во ое оо ә ф э э ә э = ото ое оо э ® э е о ов че е. 


0.4 


0.3 


E LEGEND 
| 5 F,ROLL ANGLE-RAD 





0.0 


0.0 04 0.8 1.2 1.6 2.0 2.4 2.8 
TIME-SEC 


Figure 4.20 Roll Angle vs. Time; Uncoupied Modified Roll Channel with 8, 
Centered; Continuous Open Loop System; 
Step Input: 9cz1 rad and 5645 6206509640508 86 2 0. 
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Figure 4.21 Roll Angular Rate vs. Time; Uncoupled Modified Roll Channel 
with &, Centered; Continuous Open Loop System; 
Step Input: $cz1 rad and 615620630540503 86 2 0. 


159 


LEGEND 
2 DR,ROLL TAIL INCIDENCE- DEG. 


DR,ROLL TAIL INCIDENCE- DEG. 





O 


„зз еее зд Макал ш о 


-0.4-0.3 
7 
Қ 


0.0 0.4 0.8 1.2 1.6 2.0 হি 2.8 
TIME-SEC 


Figure 4.22 Roll Tail Incidence vs. Time; Uncoupled Modified Roll Channel 
with &, Centered; Continuous Open Loop System; 
Step Input: Pe=1 rad 810 6106205306405038 ap = 0. 
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Figure 4.23 Control Surfaces deflections vs. Time; Uncoupled Modified 
Roll Channel with 8, centered; Continuous Open Loop System; 
Step Input: $cz1 rad and 610556205630640508 86 2 ০. 
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Figure 4.24 Roll Angle vs. Time; Uncoupled Modified Rol] Channel with 8. 


and ô> Centered; Continuous Open Loop System; 
бер Іприї: фгс=1 гад апо 6) -6263бас:0; 25: 0. 


142 


1.6 
1 


d: Xi 
- | 
Е .ה‎ 
‘ch | 
N | Ш 
e . 
- р: 
о | : E 
E cc EU ১ 2 ০০০১৩, 
со ј y 
E. m | 
BEER ee ו‎ -- е 
Ес Ж 
E 
E ו‎ НИВИ на S S ..Ш... 
dio 
2 T 
- | END 
Е EE -..-- ROLL ANGULAR RATE-RAD/SEC 
= ₪ 
Е \ 1 i 
N B: ; , 
ו‎ er eo nse ens MB... us... Q... 
TE E c 
0 | | 
N : : | 
° ; | : | 
eee нее ও 22 -- pete o oee a 
N ugue e: 
NE _ দে; | 
E | EI : ' 
° : : : | : | 
I | | 1 
0.0 0.4 0.8 1.2 1.6 2.0 2.4 2.8 


TIME-SĐEC 


Figure 4.25 Roll Angular Rate vs. Time; Uncoupled Modified Roll Channel 
with 6, and&, Centered; Continuous Open Loop System; 
Step Input: $051 180 810 6162063054250$ ap = 0. 
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Figure 4.26 Roll Tail Incidence vs. Time; Uncoupled Modified Roll Channel 
with ê; and ô> Centered; Continuous Open Loop System; 
Step Input:Pe=1 rad and 64,8953 —964-=93 Зв = 0. 
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Figure 4.2? Control Surfaces deflections vs. Time; Uncoupled Modified 
Roll Channel with 8,,62 centered; ContinuousOpen Loop System; 
Step Input: Фс=1 rad and645,624,63,,64:20; a, z 0. 
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Figure 4.26 Roll Angle vs. Time; Uncoupled Modified Roll Channel with &,, 
92 810 63 Centered; Continuous Open Loop System; 
560 1700: $051 180 8110 61১62063064530$ 89 2 ০, 
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Figure 4.29 Roll Angular Rate vs. Time; Uncoupled Modified Roll Channel 
with 61562 8170 ôz Centered; Continuous Open Loop System; 
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Figure 4.50 Roll Tail Incidence vs. Time; Uncoupled Modified Roll Channel 
with 61, 6, and 83 Centered; Continuous Open Loop System; 
Step Input: 9cz1 rad 8110 91062630640508 89 2 ০. 
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Figure 4.51 Control Surfaces deflections vs. Time; Uncoupled Modified 
Roll Channel : &,,5,53 Centered; Continuous OpenLoop System; 
Step Input: Yol rad and 6126206306420; ар = 0. 
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Figure 4.52 Comparison of the Roll Angle between the 
Unimpaired System and the Four Cases of Impairment. 
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Figure 4.33 Roll Angle ; Uncoupled Open Loop Channel: 
Unimpaired and Damaged Control Surfaces. 
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Figure 4.54 Comparison of the Roll Angular Rate between the 


Unimpaired System and the Four Cases of Impairment. 
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Figure 4.55 Roll Angular Rate; Uncoupied Open Loop Channel: 
Unimpaired and Damaged Control Surfaces. 
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Figure 4.36 Comparison of the Roll Tail Incidence between the 
Unimpaired System and the Four Cases of Impairment. 





154 


е 


1.1 
. 
. 

— 
. 
. 
t 
4 
. 
. 
. 
. 
. 
. 
, 
. 
. 
. 
. 
. 
. 
. 
` 
. 
. 
. 
. 
. 
. 
` 
e 
* 
e 
. 
. 
. 
- 
. 
. 
. 
. 
. 
. 
а 
. 
. 
. 
. 
. 
. 
4 
. 
. 
. 
. 
. 
. 
. 
` 
. 
. 
. 
. 
. 
. 
` 
* 
. 
* 
. 
. 
. 
` 
. 
` 
. 
. 
. 
. 
е 
. 
. 
. 
L] 
4 
L] 
4 
* 
. 
. 
. 
. 
. 
. 
% 
L 
. 
. 
. 
৪ 
. 
` 
. 
. 
. 
. 
. 
. 


1.0 


টির LEUR NUR i aca Oe sg ek ee 
сң тдейе סו‎ е Ка а 5 205 ое еее 2555৮৮5১৮১৮ =т= = ет = эе = = е = е е = р == = # Ш> ә =Ф===®а = э ёё е ө е ө ә э © ә = э ә = = = = ө ө 
1 1 


д ое еее ооо ооо оо оо е со ое мо ө ө э ое о Ё э ө дове Ќе е ৬ ক ০৩০৩০ ৬ ০০ ৩ ৯০৩৩৩০৩৬০৪৩ ৮৩০৩ ঠ৩ е ое обе ® ® е = э'ъ = = ово 


э е э е е фе o o o o o ele o 9 o ое o o» v o o o o o бо е е е до ое ооо до = = e o = одо a ® = 9 e =» в ® Фф ® ৬ э е 9 =» o o э = ao Ф = o ооо о * e 0 9 o*€w o o o 9 o o €^ a ^ 5 o > = 


ооо цесс осо ео о ৬ ৩ ও ব্রাশ ও ও ও ও ৩ ও চি ও ও গজ ৩ শি ও ৩ ৩ ০ আও ও ৬ ৩ ৩ ৩৫০ вр оо сод оо до оо цо оо е по одо то জঞ্ি 












LEGEND 
O ORIGINAL SYSTEM 
o INOPERATIVE D 1 
А INOPERATIIVE D 1,2 

+ ІМОРЕКАТІУЕ В 1,2,3 


. . . 
ооо ооо ооо ооо ооо о ও е = ৬ ও > 


ооо о ео ооо оо орао ооо оо об оо о о ә = 








0.3 0.4 0.5 0.6 0.7 0.8 9 


рото ое э э = =” =|Ф Ф > = ® = э ө э е. 


. . . 
ооо а ве ооо оса оао во аша ә ака 4 /- (5-6 өчө ә ә s э ә: эө ә өөө э ә ә өтә чө Ф эө ше эе э элее еее а ате ә ее = 
1 


-01 0.0 01 $2 


- Ж 


0.0 0.4 0.8 E 1.6 2.0 2.4 2.8 
TIME-SEC 


DR,ROLL TAIL INCIDENCE- DEG. 
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Unimpaired and Damaged Control Surfaces. 
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Y. COUPLED LINEAR PITCH AND ROLL CHANNELS 


A. INTRODUCTION 

Separate longitudinal and lateral autopilots were designed in Chapters 
Ill and IV, for the Modified Linear Open Loop system. 

In this Chapter the longitudinal and roll autopilots were combined, 
resulting in a new autopilot capable of commanding pitch, roll or combined 
maneuvers. 

The time responses to an input vector are identical to that obtained 
when the same input is put into the separate pitch and roll autopilots of the 
unimpaired system. . 

The response plots of this system are then used as a basis to 
investigate the effect of the damaged control surfaces on the combined 
system. In this section the hard over type of damage for the combined 
system is investigated. Corrective actions are then developed and presented 


and the robustness of the complete system is finally analyzed. 


B. DESIGN APPROACH, ANALYSIS AND PERFORMANSE OF THE SYSTEM 

For the design of the combined system of the circular airframe 
configuration, the requirements and assumptions mentioned in the design of 
the separate autopilots are used. To achieve the desired maneuver plane 
acceleration response for the unimpaired system, both pitch and roll 
uncoupled autopilots were designed to have the desired maneuver plane 


speed of response. In paticular, the goal for maneuver plane acceleration 


156 


response was a 0.5 second time constant for the flight condition of interest 
(j.e., 60 Kft altitude, Mach 3.95. In the next phase of the linear design 
technique the longitudinal and lateral autopilot channels are coupled using 
the aerodynamic model of Chapter III. Both channels are cross-coupled 
through the contro! surfaces (i.e., a command to the control surface of one 
channel will affect the parameters of the other and vice versa). The pitch 
and roll state Equations (III.B.5-5) through (11.8.5-17) апа (У.С.4-4) 
through (IV.C.4-13) are combined to yield a twenty-three order system 
which is presented in Appendix C. 

The Open Loop dynamic matrix is a 23 x 23 matrix. This composite 
system is formed by two block diagonal matrices as is stated in Equation 
(V.B- 1): 


| = (15א15‎ | 0 | 
ШЕ -----= L. (V.B.1-1) 
0 |Ғ2(0к10) 


That is, the upper block composed of the 13x13 F matrix from the 
pitch channel and the lower block composed of the 10x10 F matrix from the 
roll channel. 

The control distribution matrix 6, folowing the same pattern, is given 


by 8 25x10 matrix, of the form and order: 


| 6103%5) | 0 | 
6 = | -------|------ (У.В.1-2) 
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Finally the state vector X and the control input vector U matrices are 
of the form and order: 


X=- [23x1] (V.B.1-3) 


U= [10x1] (У.В.1-4) 


Because of this form, any input to the stated combined system, has a 
time response identical to that resulting from the same input into the pitch 
and roll channels separately. An example of this is shown in Figure 5.1 
where the same input vector used for the vertical translation maneuver 
longitudinally and the rolling translation maneuver laterally are combined 
and serve as inputs to the combined system resulting in a response 
identical to those achieved separately in Figures 3.9 through 3.12 and 4.12 
through 4.15. 

The combined autopilot designed in this thesis and the simulations 
made to test the different step input commands are accomplished using the 
interactive, multivariable design program called OPTSYS. 

In addition, in order to compare different control surface impairment 
configurations a computer program was developed (Appendix D). This 
program creates data files to be used as input files at the EASYPLOT 
program by changing the format of the OPTPLOT data into a fixed format 
data file. This Format is 80 characters in length with space between the 
values. 

Three different command inputs are examined for the Modified 


Unimpaired system and are discussed below. 
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Test |; For this test, both a pitch normal acceleration and a roll angle 
were commanded (1.е., М2с = 1 дее, Фе = 1 гад ап@ 61 620 630 64с = О гад). 
In Figure 5.1 the time response plots of the normal acceleration and the roll 
angle are presented. As was expected, both curves are identical to the 
corresponding plots obtained when separate autopilots were used. 

Figures 5.2 through 5.4 illustrate the performance of the other state 
variables. Figure 5.4 shows the control surface deflections for only two 
control surfaces from each channel. The remaining four coincide or are the 
mirror image of the obtained results. In the same Figure the configuration of 
the control surfaces for each channel is provided and the distinction is also 
provided in the legend. 

Test 2: For this test, only a pitch command in acceleration was applied 
(i.e, Nz, = 1 gee), while, $, = 0 rad and 845 625 635, 64, - O rad. Figures 
9.9 through 5.8 show the performance of the system. Observe that no roll 
angle or any roll rate is present, as was expected, while a downwards 
maneuver is performed according to the command input vector ( downward 
maneuver is the positive one). Figure 5.6 shows the response of the control 
surfaces for both channels. The control surfaces of the longitudinal channel 
are deflected at a +0.025 radian angle and of the lateral channel remained at 


the initial Zero position while only a downward maneuver was commanded. 
Test 3: In the last test, only a roll command is applied ( i.e., 9. - 1 rad 





), while Nz, = O gee and 845 655 630 бас = O rad. Figures 5.9 through 5.12 


illustrate the effectiveness of input vector on the state variables of 
interest. In Figure 5.12 only the two control surfaces 5, and ô, of the roll 


channel are plotted while the corresponding surfaces in the pitch channel 
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remain at zero. The other four coincide with the obtained results. The 
control surfaces of the roll channel are deflected at a +0.022 radian angle 
and of the longitudinal channel remained at the initial zero position while 
only a rolling maneuver was commanded. 

It should be mentioned here that according to the pole-Zero plot for the 
designed combined system, the combined continuous open loop system is 


also stable, since the s-plane poles are: 


6, = -109.724 + 2 (У.В-5) 
So = - 199.724 - ј18.9922 (У.В-6) 
S3 = -8.29043 + j0.05911 (Y.B-7) 
S4 = -8.29043 - j0.05911 (pitch angular rate) (У.В-8) 
Ss = -3.75933 + j2.51441 (У.В-9) 
Sp = -3.79933 - j2.51441 (У.В-10) 
S7 = -0.14393 + j0 (У.В-11) 
54 = -0.00005 + 0 (pitch tail incidence) (V.B-12) 
58 = -0.07460 + 9 (У.В-13) 
S197 -0.07460 - j0.66599 (pitch normal acceleration) (V.B- 14) 
S41 = -0.14300 + jo | (V.B-15) 
912 = -0.14300 +0 (У.В-16) 
біз < -0.14300 + 0 (У.В-17) 
$14 = -0.01850 + 10 (roll angular rate) (V.B-18) 
$15 = -9.27462 + ]29.2679 (roll angle) (У.В- 19) 
516 = -9.27462 - ]29.2679 (У.В-20) 
$17 = -2.35821 + 2.70296 (У.В-21) 
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Figure 5.1 Normal Acceleration and Roll Angle vs. Time; Combined Pitch 


and Roll Channels; Continuous Open Loop System; Step Inputs: 
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Figure 5.3 
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| га; ав = 0. 
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and Roll Channels; Continuous Open Loop System; Step Inputs 


.4 Pitch and Roll Deflections vs. Time; Combined Pitch 
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Figure 5.5 Normal Acceleration and Roll Angle vs. Time; Combined Pitch 
and Rol] Channels; Continuous Open Loop System; Step Inputs: 
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Figure 5.6 Pitch and Roll Angular Rate vs. Time; Combined Pitch | 
and Roll Channels; Continuous Open Loop System; Step Inputs: | 
619 620 630 4c = O and Noc = 1 gee, Pe = O rad; аь = 0. 
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Figure 5.7 Pitch and Roll Tail Incidence vs. Time; Combined Pitch 
and Roll! Channels; Continuous Open Loop System; Step Inputs: 
615 бәс» 630 бас = 0 апа Мес = 1 gee, Фе - 0 гад; dg - 0. 
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Figure 5.8 Pitch and Roll Deflections vs. Time; Combined Pitch 
and Roll Channels; Continuous Open Loop System; Step Inputs: 
610 6209 636 бас = О апд Мәс = [ gee, $c = 0 гад; да 2 0. 
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Figure 5.9 Normal Acceleration and Roll Angle vs. Time; Combined Pitch 
and Roll Channels; Continuous Open Loop System; Step Inputs: 
81 20 зо sc = O and Nae = O gee, Pe = 1 гад; ао = 0. 


169 


<.. fe tw mm Om a. ... Ñ-. ... ו‎ cA ה םהמ‎ оос ее ео е ае 


0.9 
1.4 


ее ено сла фо ново геого вове го ноное “ЙМ. ө ө ө ө ө ө ө Ф% е Ё ө э тә э э ө е э ө ө ө ө э ө ө э ө ө ө = ө + й ө ө = = оо ө ө э ө ө ө о ео э ө = = ө = өө „% эө өө ө ө о э ө ө ө ө ө ө ө оо ө э = ө ө ө ө ө э э ө ее ө ө ө = ө ө 
. . 


0.7 


......4--%%.“-....)ф%%-...........-Х..-.-..Ж-..............................ө................-.....................-..................... 


1.0 


. . . 
“ө..... Е«--ее«ееевееее-еФеееге«ееее»-.-егебес-е %-еегееееегее«еееееееееееееечеееееееееееееееееоеееекееееееовеееоеееееовкееееееееекек-е-еееее.--ғ4.... 
: . 


0.6 


0.8 


ано Да ененнен е ноар 1-2 -...-1-..11 בש‎ DS  " 


0.6 


0.4 


вве агане ב בב‎ 222 a A 


PILCH ANGULAR RATE-RAD/SEC 


0 
P,ROLL ANGULAR RATE-RAD/SEC 
0.4 


0.3 


1 
=+ a amas 5 чоя ж жоя чоя чо су 1 =ош ө Ф =ош ө ө =ош = = ө = = ө = = ө ө = = оф ө 5 Ф ө = 5 э а + э +}ь „ „ = = & & . 5 ө э ө ө ө ө ө ө ө ө ө ө э ө ө ө ө ө = ө ө = э „фр ө тө ө ө э э ө = Өө ө ө ө э = ө ө ө *® ө © = =ф ө Өө Ө = Ө ө өө э э Фо өө ө ө Ө ө э о 





оз 
a © LEGEND 
e 
e : 
T Ө রা রি EE e = - эй 
2: ইউ. E LTEM MAUS T 
—_ — 
а te 
> са 
о a | ы M 


0.0 a 0.8 1.2 6 . 2.0 | с 
TIME - SEC 


Figure 5.10 Pitch and Roll Angular Rate vs. Time; Combined Pitch 
and Roll Channels; Continuous Open Loop System; Step Inputs: 
6ו6‎ 520 530 546 = 0 and Noe = 0 gee, Pe = 1 rad; ae = 0. 
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Figure 5.11 Pitch and Roll Tail Incidence vs. Time; Combined Pitch 
and Roll Channels; Continuous Open Loop System; Step Inputs: 
816 620 630? бас = 0 апа Мәс - 0 дее, Фе - 1 rad; Зе = 0. 
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12 Pitch and Roll Deflections vs. Time; Combined Pitch 
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518 = -2.59821 - |2.20296 (У:В-22) 


$13 = -8.89209 + jO (৬.৪-25) 
S2g = -9.18517? + jO (roll tail incidence) (V.B-24) 
$21 = -15.0000 + ]0 (V.B-25) 
522 = -15.0000 + 0 (V.B-26) 
$95 = -15.0000 + jo (V.B-27) 


C. THEEFFECTS DUE TO ACTUATOR IMPAIRMENT 

In this section the effects on the system due to a damaged actuator 
was investigated. One control surface (6,) was considered damaged. It was 
locked at a fixed angle while the other three control surfaces had a variable 
deflection angle. 

Two different control inputs were examined for this type of 
impairment and are discussed below. For each test, corrective actions were 
then developed and presented. 

Test 1. Command Input: 

(1) 6 for both channels locked at 0.2 rad 





(2) 55, 53, 54 for both channels have variable deflections 
(3) Мес 1 дее 
(4) Фс = 0 гай 
Using the above control input the time response plots of the system 
were obtained, analyzed and presented in Figures 5.13 through 5.16. In 


Figure 5.13 the pitch normal acceleration and the roll angle are presented. 
Commanding a downward maneuver (i.e., Nz =1 gee and $ = 0) the missile 
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body the first tenth of the second moves downward and then upward with its 
pitch axis directed at the velocity vector until it reaches the desired 
maneuver level. This is due to the induced pitching moment effect of the 
locked 8, control surface and the undamaged control surfaces (85, 85, 64} 
which try to counteract the impairment of 8,. Allthough no roll motion is 
required to maintain cordination for this maneuver, as is shown in Figure 
5.13 а roll angle (0.4 rad) is present. In Figure 5.14 the pitch and roll 
angular rates are also presented in which the 1.9 rad/sec roll rate is 
achieved. In Figure 5.16 two control control surfaces from each channel are 
presented. The other two are the mirror images of the obtained results. The 
ô, (of pitch and roll channels) respond properly to the command input. The 
ô> for the pitch channel is deflected 0.174 radians and for the roll channel 
0.130 radians. Both are stabilized after 1.4 seconds. 

Corrective actions: The roll effect of the damaged control surface on 
the performance of the combined autopilot was examined here and 
corrective actions to overcome this problem were studied and presented. 

Figure 5.13 shows the roll angle (9) reaching its final value of 0.40 rad 
(239) in approximately 1.4 seconds with a 7% overshoot. The first try to 
overcome the roll effect was to command the opposite surface (85) in an 
equal angle of deflection 0.2 radian. The obtained time response of the pitch 
normal acceleration and the roll angle are presented in Figure 5.17. For 
comparison, Figure 5.17 presents three configurations. It shows that the 
produced roll angle due to the 8, control surface impairment is eliminated. 
But, there is a penalty in the pitch normal acceleration which reaches 


approximately 4 gee and a oscillatory behavior is evident before the steady 
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state value of 1 gee command being achieved. So, another test was 
performed by commanding the &3 surface at an opposite angle of deflection 
(-0.2 rad). By this action, shown in Figure 5.18, the normal acceleration 
responds very well but the roll angle is increased at 1.18 radians (689) after 
1.75 seconds, as might be expected. The next test was performed using the 
second control surface (85). In Figure 5.19 the time response of the same 
variables of the system are obtained and presented. The control surface ë; 
is commanded at the same angle of deflection i.e., 0.2 radians. As is 
illusrtated in Figure 5.19 commanding a downward maneuver the missile 
body reaches the desired maneuver level without any roll motion (i.e., ?א‎ =[ 
дее апа % = 0). 

At the end of this Chapter a transition reconfiguration strategy is 
proposed addressing the topics of failure detection and isolation, and 
implementation. That is, one. pssible strategy for transitioning from the 
normal control laws for a missile with all control surfaces operative, to 
one set of reconfigured control laws for a missile with a failed surface. 

Test 2. Command Input: 

(1) 8, for both channels locked at 0.2 rad. 

(2) 85, 63, 64 for both channels have variable deflections. 
(3) Nz, = 0 gee. 

(4) $, = 1 rad. 

Using the above step input vector the time response plots of the sustem 
were obtained and presented in Figures 5.20 through 5.25. 

In Figure 5.20 the time response of the pitch normal acceleration and 


the roll angle of the combined system are illustrated. Due to the fact that 
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the two channel autopilots were cross-coupled through the four control 
surfaces (81, 62, 63, 64), a damage in any actuator(s) will affect the 
performance of both channels. The locked control surface 6, at the one rad 
angle affects the combined autopilot. A constant incorrect level was 
achived іп the response of the roll angle (i.&, 9 = 1.4 rad instead of the 
commanded one rad). In addition, due to the induced pitching moment effect 
by the damaged control surface, a down and upward maneuver was 
commanded bu the pitch channel during the first second of the response of 
the normal acceleration. 

A change of the other state variables was obvious as is shown in 
Figures 5.21 through 5.23. Figure 5.25, which represents the time response 
of two control surfaces from each channel, i.e., 8, and 85, indicates the 
commanded (simulated) failure of 0.2 radian at 6, and the deflection angle 
of 0.13 radians (7.59) for the other three control surfaces. The control 
surfaces not presented in this plot are the mirror images of the obtained 
ones. 

Corrective actions: The roll and pitch effects caused by the damaged 
control surface ô; on the combined autopilot performance was examined and 
corrective actions to overcome this failure were studied and presented. 

Іп Figure 5.20 the undesirable rolling and pitching effect caused бу 
the damaged 8, control surface is obvious. The missile rolls about its 
velocity vector with a 1.4 radians roll angle (809) (instead of the 
commanded 1 radian). Also observe the effect on the pitch normal 
acceleration which should be constant at zero since no down or upward 


maneuver was Commanded. 
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To overcome these failures several runs of the Controls Program 
OPTSYS were performed to find an appropriate corrective action for the 
control system. The most interesting cases are illustrated in Figures 5.24 
and 5.25. 

In Figure 5.24 the roll angle (9) and the pitch normal acceleration (n;) 
were obtained and presented for three cases to compare their responses. À 
comparison was performed between the original system, the impaired 
system with the ô; effector locked at an angle of 0.2 radians, and the 
system with the ô, control surface commanded at an angle equal but 
opposite (-0.2) to the defective surface. As is illustrated in Figure 5.24, if 
such an input was commanded, the missile body rolls even more at a 2.2 
radians (126°) angle as might be expected. Also, the pitch normal 
acceleration increases and achieves the large value of -4.2 gee's before 
stabilizing after 1.4 seconds. И 

In Figure 5.25 the roll angle (9) and the pitch normal acceleration (n;) 
were obtained and presented but now the second control surface 52 was 
commanded to be deflected at the same angle as the §,, (0.2 radians). As is 
illustrated in Figure 5.25, both effects, on the roll angle and the pitch 
normal acceleration, were eliminated and their time response plots 
indicates that the missile body performed according to the command input. 

In comparing output time responses, as was done by these two tests, it 
is apparent that the impaired missile was able to track the unimpaired 
missile as long as there was another control surface to provide the 


appropriate corrective action and control power of the failed surface. 
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In the case where one or more control surfaces are inoperative or shot 
away and are not considered here, similar effects are expected. Thís is due 
to the fact that as far as all the control surfaces .are operative, are 
commanded by the system such that one will counteract the induced moment 
effect of the other. But when one control surface is inoperative or shot 


away undesirable roll, pitch or combined effects will be present. 


D. PROPOSED RECONFIGURATION STRATEGY 
1. Introduction 
This section presents one possible strategy for transitioning from 
the control laws for a missile with all surfaces operative, to a set of 
reconfigured control laws for an missile with a failed surface. The 
discussion begins with a definition of reconfiguration and why it is 
important. A transition stategy is then proposed, addressing the topics of 
failure detection, isolation, and implementation. 
2. Reconfiguration 
Reconfiguration is a technique which shows promise in restoring 
missile stability and a degraded level of controllability in the event a 
control surface becomes inoperative. This technique implies that the control 
laws governing missile motion compensate for the failed surface by using 
the remaining surfaces to maintain control. The control laws of Chapter III, 
IV and the corrective actions of Chapter V provide the independent control 
of each surface necessary to implement reconfiguration which deals with 


the coupled aerodynamics. 
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3. Transition Strategy 


The development of reconfigurable control laws, such as those 
developed in this thesis, must include a method of failure identification and 
implementation of a revised control law. For a given failure it thus 
necessary to determine when the new laws should be used and how to 
switch to them from from normal control laws. A failure detection and 
isolation scheme must be developed to determine when and which surface 
has failed. The proposed failure detection and isolation system would use a 
set of redundant sensors and an estimator, to determine if the actuators are 
Still operative. Should an actuator fail to move a control surface, the 
sensors notify a control surface monitor. The monitor validates the fail 
signal and determines which control surface is inoperative. 

The control surface monitor is the heart of the reconfigurable flight 
control system. When notified of a failure, the monitor informs the flight 
control computers of the problem and provides the necessary information 
concerning which set of new control laws (corrective actions) should be 
used to calculate the appropriate surface deflections. The off-line flight 
control computer is programmed with the new control laws and then signals 
the on-line system that it is ready to assume control of the missile. Thus 
the autopilot is again ready to control a stable missile with some degree of 
control capability. 
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Normal Acceleration and Roll Angle vs. Time; Combined Pitch 
and Roll Channels; Continuous Open Loop System; Step Inputs: 
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Figure 5.14 Hard Over $, Type of Impairment: 
Pitch and Roll Angular Rate vs. Time; Combined Pitch 
and Roll Channels; Continuous Open Loop System; Step Inputs: 
&1650.21840, бос» 630 бас - 0 апа Мәс = | gee, Фе = 1 rads дь = 0. 
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Figure 5.15 Hard Over 5, Type of Impairment: 


Pitch and Roll Tail Incidence vs. Time; Combined Pitch 
and Roll Channels; Continuous Open Loop System; Step Inputs: 
61050.2180, 620 630) 640 3 0 апа №. = 1 дее, ф. = 0 габ; а, = 0. 
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Figure 5.16 Hard Over 5, Tupe of Impairment: 
Pitch and Roll Deflections vs. Time; Combined Pitch 
and Roll Channels; Continuous Open Loop System; Step Inputs: 
61c70.2rad, 620 635, 6425 0 800 No, = 1 gee, Pe = О rad; a, = 0. 
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Figure 5.1? First Corrective Action and Comparison of the Pitch 
Norma) Acceleration (Nz) and the Roll Angle (9), between 
Three Different Effectors Configuration. 
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Figure 5.18 Second Corrective Action and Comparison of the Pitch 
Normal Acceleration (Nz) and the Roll Angle (9), between 
Three Different Effectors Configuration. 
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Figure 5.19 Final Corrective Action and Comparison of the Pitch 
Normal Acceleration (Nz) and the Roll Angle (9), between 
Three Different Effectors Configuration. 
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Figure 5.20 Hard Over &, Type of Impairment: 
Normal Acceleration and Roll Angle vs. Time; Combined Pitch 
and Roll Channels; Continuous Open Loop System; Step Inputs: 
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Figure 5.21 Hard Over 86, Type of Impairment: 
Pitch and Roll Angular Rate vs. Time; Combined Pitch 
and Rall Channels; Continuous Open Loop System; Step Inputs: 
51г:20.2г20, 6-,, бе бас: 0 and No, = O gee, P, = 1 rad; a, = 0. 
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Figure 5.22 Hard Over &, Type of Impairment: 
Pitch and Roll Tail Incidence vs. Time; Combined Pitch 
and Roll Channels; Continuous Open Loop System; Step Inputs: 
5 1с=0.2гад, 5 бас, 64, = 0 and Nop = 0 gee, $e = 1 rad; a, = 0. 
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Figure 5.23 Hard Over &, Type of Impairment: 
Pitch and Roll Deflections vs. Time; Combined Pitch 
and Roll Channels; Continuous Open Loop System; Step Inputs: 
§1,=0.2rad, Son, 6305 64, 2 O and No, = 0 gee, P, = 1 rad; a, = 0. 
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Figure 5.24 First Corrective Action and Comparison of the Pitch 
Normal Acceleration (Nz) and the Roll Angle (9), between 
Three Different Effectors Configuration. 
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Figure 5.25 Second Corrective Action and Comparison of the Pitch 


Normal Acceleration (Nz) and the Roll Angle (9), between 
Three Different Effectors Configuration. 
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E. ROBUSTNESS OF THE UNIMPAIRED SYSTEM 
1. Introduction 

Up to this point the continuous open loop system was investigated 
and it was found to be stable, as is indicated by Equations (V.B-5) through 
(V.B-27). In this section the concept of feedback is introduced and the 
robustness of the system is examined. 

Design criteria for feedback control systems can be lumped into 
three broad categories: performance, stability, and robustness. 
Performance criteria are generally specified as a requirement for the 
control system to follow commands or reject disturbances. The most 
commonly used stability criteria are stability margins which generally 
refer to the amount of gain or phase variations from a design condition a 
system employing feedback control can experience before becoming 
unstable. 

in early feedback control applications, reducing sensitivity was 
the only control design criteria, particularly reducing sensitivity to noise. 
Over the last several years, sensitivity issues have attracted new attention 
and have generated a new requirement, robustness. A control system is 
robust (insensitive) if performance is maintained in the presence of 
variations in the plant from the design model. These variations can be 
modeling approximations, modeling uncertainties, failures of control 
elements (sensors, actuators), noise, non-linearities, etc. 

These criteria are not mutually exclusive. Good stability 


characteristics in well designed systems are usually accompanied by good 
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performance. Robustness must be considered relative to the variations from 
the design model that the control system will! experience. 


2. The Benefits of Feedback and the Singular Value Desiqn Technique 








Figure 5.26 Standard Feedback Configuration 


We deal with the standard feedback configuration illustrated in 
Figure 5.26. It consists of the interconnected plant (G) and controller (K) 
forced by commands (r), measurement noise (n), and disturbances (d). 

The noise in this case is assumed to be zero. The disturbances are 
assumed to be reflected to the measured output (y). 

All signals are multivariable, in general, and both nominal 
mathematical models for G and K are finite dimensionali linear time 
invariant systems with transfer functions matrices G(s) and K(s). Then it is 
well known that the configuration, if it is stable, has the following major 
properties: 


а Input-output behavior 


у=6К (+ 6К)-1 (г -п) + (+6К)-1 9 (V.E.2-1) 
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e=r-U 
zOI*«G6K-!(ir-d*«GK(I*GK)!n (V.E.2-2) 


b. System sensitivity 
АН а= (1+ 6' К)-! АН у (У.Е22-3) 


In Equation (V.E.Z-3), AH „ү and AH gi denote changes in the closed 
loop system and changes in a nominally equivalent open loop system, 
respectively, caused by changes in the plant 6, i.e., 6' 2 G * AG. 

Equations (V.E.2-1) through (V.E.2-3) summarize the fundamental 
benefits and design objectives inherent in feedback loops. Specifically, 
Equation (V.E.2-2) shows that the loop's errors in the presence of commands 
and disturbances can be made "small" by making the sensitivity operator, 
(I + G KYI, "small". Equation (V.E2-3) shows that loop sensitivity is 
improved under these same conditions, provided AG is not too large 
from G. 

With the increased interest in MIMO systems numerous methods of 
design һауе been employed to obtain suitable system performance and 
robustness with varying degrees of success. One primary method of design 
is to keep the plant as decoupled as possible throughout the design so that 
each individual element may be controlled independently and designed 
essentially as a single loop system. In another procedure the multiloop 
system is modified into a system that has diagonal elements that are much 


larger than any off-diagonal elements. This diagonally dominant system is 
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then in a form where conventional Nyquist type techniques can be employed 
in the analysis. A third common MIMO design method is that of the Linear 
Quadratic (LQ) method. This method uses a quadratic cost functional and 
optimization principles to allow the designer to design for various 
performance levels by adjusting the matrix weighting terms used in the 
cost function. The major difficulty with all of the above methods is that 
they are not necessarily robust. This is especially true for cross-coupling 
terms between loops. 

As was mentioned before, a simple interpretation of robustness is 
the ability of the system to tolerate design pertubations. These 
perturbations could be in the form of actuator failures, plant parameter 
uncertainty, unmodeled dynamic or nonlinear terms, or any one of many 
other perturbations to the nominal design of the system. 

The analysis of the robustness of the system was obtained using 
the Pole Placement and Robustness Design Program (POPLAR) developed by 
Gordon, V.C. [Ref. 5]. 

Utilizing the above program and using the data file shown in 
Appendix E for the state-feedback autopilot, the unimpaired system was 
found to be robust. The standard full state-feedback design was carried out 
using the N.P.S. program OPTSYS. 

Figures 5.27 and 5.28 show the Minimum Additive Output Singular 
value (SVADMO) vs. frequency, and the Minimum Additive Input Singular 
Value (MIN.ADD.IN.SV) vs. frequency respectively for the state-feedback 
designed autopilot. It is noted that for 0.1 to 100 db frequency the output 
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F*®G singular value is too low ( 0.03 rad ) but the input transfer function 
G*F input minimum singular value is above the 0.5 rad region. 

Therefore the closed loop unimpaired system is robust (at least in 
the input). Further analysis would be needed if robustness were required in 
both the input and output case simultaneously. 

In the case of state-feedback autopilot and for a singular value of 
0.5 rad the pole placement and robustness design routine places the poles as 
follows: | 


Ordered (input) Computed (output) 


-1706.00000 0 -1706.58980 0 
-285.39990 0 -285.42944 0 
-219.48999 0 -158.78516 0 
-117770000 0 -49.74738 0 

-56.00000 0 -49.74738 0 

-21.84399 0 -27.69525 0 
-7.05800 -j1.59000 -27.59982 0 
-7.05800 0 -8.57168 -|2.58016 
-4970000 0 -857166 6 
-4970000 0 -6.42040 0 
-4.88000- 0 -6.35900-|5.27725 
-4880000 0 -6.35900 5 
-2.90000 0 -5.99561 0 
-0.77000 0 0 -427422  j0.00000 
-0.47000 0 -2.33516 0 
-0.47000 0 -098429 0 
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Figure 5.27 SVADMO vs. Frequency; Copled Pitch and Roll Channel 
Autopilot; State-Feedback Design; Continuous Closed 
Loop Unimpaired System; Circular Airfram. 
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Figure 5.28 MIN ADD IN SY vs. Frequency; Copled Pitch and Roll Channel 
Autopilot; State-Feedback Design; Continuous Closed 
Loop Unimpaired System; Circular Airfram. 
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VI. CONCLUSIONS AND RECOMMENDATIONS 


CONCLUSIONS 


The purpose of this thesis was to examine the control deficiecies in a 


Bank-to-Turn (B.T.T) cruise missile, due to specific impairments at its 


control surfaces during the intermediate phase of its mission. The following 


are the principal conclusions based on this work. 


(1) The original and modified classical designed autopilots which involve 


separate inputs for each control surface were proved to have 
identical performance for the longitudinal and lateral channels. 


(2) The investigation of specific type of damages at the B.T.T. missile's 


control surfaces, showed that the hard over case ( i.e., one control 
surface locked at a certain deflection angle) shows greater excursion 
in the normal acceleration and roll angle with large values being 
achieved, than the other type of impairments. 


It was showed that the impaired missile was able to track the 
unimpaired missile as long as there was another control surface to 
provide the appropriate corrective action and the control power of the 
failed surface. 


A discussion was given of a possible method for transitioning from 
normal control laws to a set of reconfigured control laws for a 
missile with failed surface. The strategy requires the system to 
detect and isolate a surface failure and then transfer control to an 
off-line flight control computer which has been programmed with the 
new control laws (corrective actions). 


(3) 


(4) 


(5) The state-feedback autopilot was introduced as additional dynamic 


design in order to implement control to the unimpaired system. This 
design, analyzed in terms of its transient response, was found to 
meet the desired requirements. 
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(6) The performance of the coupled pitch and roll channel was found to be 
satisfactory and the unimpaired system proved to be robust. 


B. RECOMMENDATIONS 


Finally the following recommendations should be taken into 
consideration: 


(1) More complex and other type of impairments could be investigated 
using the designed autopilot, in order to examine the overall 
controllability of the missile. 


(2) In order to make the system simpler, some of the returning gain loops 
of the state-feedback design must be eliminated. A further 
investigation must then be conducted in order to examine if the 
performance of the resulting design remains robust. 


(3) Further analysis should be conducted if robustness were required in 
both the input and output case simultaneously. 
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APPENDIX A 


PLANT SYSTEM AND INPUT MATRICES OF CONTINUOUS OPEN LOOP 
PITCH CHANNEL AUTOPILOT (UNIMPAIRED SYSTEM) 


OPEN LOOP DYNAMICS MATR1X................. Е 
-0.15000+03 — 0.0000D«00  ০*০০০০০+০০ | 0.00000*06 0.00000900 0.00000%00 0.0000D+00 0.00000*00 0.00000+00 0.00000%00 
0.00000%00 0.00000»00 3 Е 
9.1353D*01 0.00000•00 0.10000°01 0.00000*00 0.0000000 0.00000*00 0.0000D+00 0.00000%00 0.00000+00 0. 0 
0.00000+00 0.00000%00 0.00000%00 
-0.65720*01 0.00000%00 -0.5000001 0.00000%00 0.00000%00 0.00000%00 0.00000%00 0.00000+00 0.00000%00 0.00000%00 
0.00000ғ00 0.0000000 0 
0.0000D*00 0.00000%00 0.0000D+00 0.0000000 | 0.10000+01 0.00000+00 0.00000%00 -0.59260%02 0.00000%00 0.00000%00 
0.00000%00 0.00000+00 0.00000+00 
0.00000%00 0.00000+00 0.00000+00 -0.4436D*02 -0.14920+00 0.00000%00 0.0000000 0.17890%01 0.00000+00 0.00000%00 
0.0000000 0.00000+00 0.00000%00 
0.00000%00 0.00000%00 0.00000%00 0.00000»00 0.0000D*00 0.0000000 0.10000%01 -0.70190%00 0.00000%00 0.00000%00 
0.00000%00 0.00000%00 0.0000D+00 
0.00000%00 0.00000+00 0.00000%00 0.0000000 0.00000%00 -0.6436D*02 -0.1492D«00 -0.1043D*04 0.00000%00 0.00000%00 
0.0000D*00 0.00000+00 0 
0.00000%00 0.00000»00 0.00000%00 0.0000D*00 0.0000D*00 | 0.00000+00 | 0.00000+00 -0.18840%05 -0.94200%02 0.96200%02 
-0.94200+02  0.9420D*02 0.00000%00 | 
0.2306D*00 0.11150%01 0.1703D*00 -0.1115D*01 -0.17030+00 0.00000%00 0.0000000 0.10090%02 -0.14500%00 0.00000%00 
0.0000000 0.0000000 0.00000%00 
-0.2304D*00 -0.1115D*01 -0.1703D*00  0.1115D*01 0.1703D*00 0.0000D*00  0.0000D*00 -0.10090%02 0.00000+00 -0.14300*00 
0.00020%00 0.00000+00 0 
0.2304D*00 | 0.11150+01 0.17030+00 -0.11150+01 -0.17030%00 0.00000%00 0.0000D+00 0.10090%02 0.00000%00 0.00000%00 
-0.1420000 0.00000+00 0 
-0.2304D*00 -0.11150%01 -0.17050%00 0.11150%01 0.17050%00 0.00000%00 0.00000+00 -0.10090+02 0.00000%00 0 
0.00000+00 -0.14500%00 0.00000+00 
0.0000D*00 0.0000D*00 0.0000000 0.0000000 0.00000%00 0.00000%00 -0.10000901 0.88710+05 0.46320%05 -0.66320%02 
0.44520+0$ -0.46432D*03 0.00000%00 
ОРЕМ (ООР EIGENVALUES.......... DETCSI-F).. 
Ш. 597 + ЊЕ 
26E*02, 1.89922E+01: 8.29043E+00, 8.05911E+00:-3.75933E+00, 2.51441E+00: 


1 .-..+ 


1.43000Е-01. 2 


6.65991Е+00:-1.43000Е-01:-1.45000Е-0] : 


THE CONTROL DISTRIBUTION MATRIX.......... G... 
0.0000D+00 0.0000D+00 0.0000D+00 0.0000D+00 0.0000D+00 
0.0000D+00 0.0000D+00 0.0000D+00 0.00000+00 -0.1353D+01 
0.0000D+00 0.0000D+00 0.0000D+00 0.0000D+00 0.6572D+01 
0.0000D+00 0.0000D+00 0.0000D+00 0.0000D+00 0.0000D+00 
0.0000D+00 0.0000D+00 0.0000D+00 0.0000D+00 0.0000D+00 
0.0000D+00 0.0000D+00 0.0000D+00 0.0000D+00 0.0000D+00 
0.0000D+00 0.0000D+00 0.0000D+00 0.0000D+00 0.0000D+00 
-0.9420D+02 0.9420D+02 -0.9420D+02 0.9420D+02 0.0000D+00 
-0.1430D+00 0.0000D+00 0.0000D+00 0.0000D+00 -0.2304D+00 
0.0000D+00 -0.1430D+00 0.0000D+00 0.0000D+00 0.2304D+00 
0.0000D+00 0.0000D+00 -0.1430D+00 0.0000D+00 -0.2304D+00 
0.0000D+00 0.0000D+00 0.0000D+00 -0.1430D+00 0.2304D+00 
0.4432D+03 -0.4432D+03 0.4432D+03 -0.4432D+03 0.0000D+00 
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CONTROL MATRIX...C=BINVxGT%xS... 


THE OPTIMAL FEEDBACK GATN 


1.45135р-01. -8.1227р=01 =1 54932070F 1.57060-02 5.235622-05 -2. 3920 В p 
9.31210202 -9 5242900 -5 18 БЕТІ 
-1.4313D-01 8.1227D-01 1.56522-0[ 7-1557064D-00 925. 2562D- 0 2 .1 
৮9772761725 Оо.51210=02 8 .1 
1.6313530-01 -8.1227D-01 -1.5452D-01 1.57040-02 5.2562р-053 -2.4426р-01 
9.51212-02  -9.51210502 29 1 
-1.6515D=01 8.1227D-01 1.5652D-01 -1.5704D-02 -3.2362D-03 2.6626D-01 
৮5021915111 9esT21D-02 8.1761D-01 
-5.5195-03 -2.566%Вт-00---46,55З5 Рр ПР Тора 6 .67 2 1.78600-01 
-3.89500-01 5.89500-01 9.2606D-02 


C-LOOP OPTIMAL REG. E-VALUES...DETCSI-F+GxC).. 


-9.09152Е%02:-2.11195Е%02:-1.78501Е%01:-9.7953506Е%00:-5.22654Е%00:-6.87901Е-01: 
ОЕ 6.65991E+00:-1.46233E-01:-5.42351E-06:-1.43000E-01:-1.43000E-01: 


₪ c s, T———— НИ 
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APPENDIX B 


PLANT SYSTEM AND INPUT MATRICES OF CONTINUOUS OPEN LOOP 
ROLL CHANNEL AUTOPILOT (UNIMPAIRED SYSTEM) 


THE SYSTEM MATRIX "F"-MATRIX ... ` 
0.00000 0.00000 . 0.00000 0.00000 0.00000 0.00000 
0.00000 0.00000 0.00000 8.05662 
1.00000 0.00000 0.00000 0.00000 0.00000 0.00000 
0.00000 0.00000 0.00000 0.00000 
0.00000 -17.60000 -8.00000 0.00000 0.00000 0.00000 
0.00000 0.00000 0.00000 0.00000 
80. 75525 “088550 5 -5.00000 0.00000 0.00000 
0.00000 0.00000 0.00000 -0.40554 
0.00000 0.00000 0.00000 0.00000 -6.00000 0.00009 
0.00000 0.00000 0.00000 0.68291 
29710271 -0.12048 0.04795 14.31820 -14.18184 -15.00000 
0.00000 0.00000 0.00000 -0.14067 
0810271 0.12048 -0.04795 -14.31820 14.18184 0.00000 
-15.00000 0.00000 0.00000 0.14067 
0.10271 0.12048 -0.04795 -14.31820 14.18184 0.00000 
0.00000 -15.00000 0.00000 0.14067 
=Ш 10271 -0.12068 0.04795 14.31820 -19.18184 0.00000 
0.00000 0.00000 -15.00000 -0.14067 
0.00000 0.00000 0.00000 0.00000 0.00000 2698.83000 
-2698.85000 -2698.85000 2698.83000 -188.40000 
EI LOOP EIGENVALUES.......... ב‎ 


~1.85057E+02:-9.27462E+00, 2.92679E+01:-2. 35821E+00, 2.70296Е%00:-8 


; + 
-5.18517E*00:-1.50000E*01:-1.50000E*01:-1.50000E*01: со Ес 
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THE OPTIMAL 
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APPENDIX C 


PLANT SYSTEM AND INPUT MATRICES OF CONTINUOUS OPEN LOOP 
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0.00000+ 00 
0.0000D*00 


0.0000D*00 
0.0000D*00 


0.0000D*00 
0.0000D*00 


0.0000D*00 
0.0000D*00 


0.0000D*00 
0.0000D*00 


0.0000D*00 
0.00000*00 


0.94200%02 
0.00000%00 


0.00000%00 
0.00000%00 


-0.14300%00 
0.00000%00 


0.00000%00 
0.00000%00 


0.0000D*00 
0.0000D*00 


-0.6632D*03 
0.0000D*00 


0.0000D*00 
0.0000D*00 


0.0000D*00 
0.0000D*00 


0.0000D*00 
0.0000D*00 


0.0000D*00 
0.0000D*00 


0.0000D*00 
0.0000D*00 


0.0000D*00 
0.0000D*00 


0.0000D*00 
-0.1500D*02 


0.0000D*00 
0.0000D*00 


0.0000D*00 
0.0000D*00 


0.0000D*00 
-0.2699D*04 


0.0000D*00 
0.0000D*00 


0.0000D*00 
0.0000D*00 


0.0000D+00 
0.00000%00 


0.00000%00 
0.00000%00 


0.0000D*00 
0.0000D*00 


0.0000D*00 
0.0000D*00 


0.0000D*00 
0.0000D*00 


-0.94200%02 
0.0000D*00 


-0.14300%00 
0.0000D*00 


0.0000D*00 
0.0000D*00 


0.0000D*00 
0.0000D*00 


0.0000D*00 
0.0000D*00 


0.44 520203 
0:00000%00 


0.00000%00 
0.0000D*00 


0.0000D*00 
0.0000D*00 


0.0000D*00 
0.0000D*00 


0.0000D*00 
0.0000D*00 


0.0000D*00 
0.0000D*00 


0.0000D*00 
-0.1500D*02 


0.0000D*00 
0.0000D*00 


0.0000D*00 
0.0000D*00 


0.0000D*00 
0.0000D*00 


0.0000D*00 
0.2699D*06 


0.0000D*00 
0.0000D*00 


0.0000D*00 
0.0000D*00 


0.0000D*00 
0.0000D*00 


-0.59260%02 
0.0000D*00 


0.1789D*01 
0.0000D*00 


-0.7019D*00 
0.0000D*00 


-0.1063D*06 
0.0000D*00 


-0.18840%03 
0.0000D*00 


0.1009D*02 
0.0000D*00 


-0.1009D*02 
0.0000D*00 


0.1009D*02 
0.0000D*00 


-0.1009D*02 
0.0000D*00 


0.8871D*03 
0.0000D*00 


0.0000D*00 
0.0000D*00 


0.0000D*00 
0.0000D*00 


0.0000D*00 
0.0000D*00 


0.0000D*00 
0.0000D*00 


0.0000D*00 
-0.6000D*01 


0.0000D*00 
-0. 14180+02 


0.0000D*00 
0.6182 


0 0 
0.2 


0.00000%00 
-0.14180%02 


0.0000D*00 
0.00000»00 


0.00000%00 
0.0000D*00 


0.0000D*00 
0.0000D*00 


0.0000D*00 
0.0000D*00 


0.0000D*00 
0.0000D*00 


0.0000D*00 
0.0000D*00 


0.1000D*01 
0.0000D*00 


-0.1692D*00 
0.0000D*00 


0.0000D*00 
0.0000D*00 


2.00000%00 
0.0000D*00 


0.0000D*00 
0.0000D*00 


0.0000D*00 
0.0000D*00 


0.0000D*00 
0.0000D*00 


-0.1000D*01 
0.0000D*00 


0.0000D*00 
0.0000D*00 


0.0000D*00 
0.0000D*00 


0.0000D*00 
0.0000D*00 


0.0000D*00 
-0.5000D*01 


0.0000D*00 
0.0000D*00 


0.0000D*00 
0.14320%02 


0.0000D*00 


-0.14320%02 


0.00000%00 


-0.14520%02 


0.00000%00 
0.14320+02 


0 0 


0.0000D*00 


0.0000D*00 
0.0000D*00 


0.0000D*00 
0.0000D*00 


0.0000D*00 
0.0000D*00 


0.0000D*00 
0.0000D*00 


0.0000D*00 
0.0000D*00 


0.0000D*00 
0.0000D*00 


-0.6636D*02 
0.0000D*00 


0.0000D*00 
0.0000D*00 


0.0000D*00 
0.0000D*00 


0.0000D*00 
0.0000D*00 


0.0000D*00 
0.0000D*00 


0.0000D*00 
0.0000D*00 


0.0000D*00 
0.0000D*00 


0.0000D*00 
0.0000D*00 


0.0000D*00 
0 0 


0.0000D*00 
-0.8000D*01 


0.00000%00 
0 


0.00000%00 
0.00000%00 


0.00000%00 
0.47950-01 


0.0000D*00 


-0.6795D-01 


0.0000D*00 


-0.47950-01 


0.00000%00 
0.47950-01 


0 0 
0.0 


0.00000%00 
0.00000%00 


0.00000%00 
0.0000D*00 


0.0000D*00 
0.0000D*00 


0.1000D*01 
0.0000D*00 


-0.14920+00 
0.00000%00 


0.0000D*00 
0.0000D*00 


0.00000%00 
0.00000%00 


0.00000%00 
0.0000D*00 


-0.1703D*00 
0.0000D*00 


0.1703D*00 
0.0000D*00 


- 90 
0.00000%00 


0.1703D*00 
0.0000D*00 


0.0000D*00 
0.0000D*00 


0.0000D*00 
0.0000D*00 


0.0000D*00 
0.0000D*00 


0.0000D*00 
-0.1760D*02 


0.0000D*00 


-0.88550+00 


0.0000D*00 
0.0000D*00 


0.0000D*00 


-0.1205D*00 


0.0000D*00 
0.1205D*00 


0.0000D*00 
0.1205D*00 


0.0000D*00 


-0.1205D*00 


0.0000D*00 
0.0000D*00 
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0.0000D*00 
0.0000D*00 


0.0000D*00 
0.0000D*00 


0.0000D*00 
0.0000D*00 


0.0000D*00 
0.0000D*00 


-0.6636D*02 
0.0000D*00 


0.0000D*00 
0.0000D*00 


0.0000D*00 
0.0000D*00 


0.0000D*00 
0.0000D*00 


- 1 
0.00000%00 


0.1115D*01 
0.0000D*00 


-0.1115D*01 
0.0000D*00 


0.1115D*01 
0.0000D*00 


0.0000D*00 
0.0000D*00 


0.0000D*00 
0.0000D*00 


0.0000D*00 
0.1000D*01 


0.0000D*00 
0.0000D*00 


0.0000D*00 
-0.75330%00 


0.0000D*00 
0.0000D*00 


0.0000D*00 
-0.1027D*00 


0.0000D*00 
0.1027D*00 


0.0000D*00 
0.1027D*00 


0.0000D*00 
-0.1027D*00 


0.0000D*00 
0.0000D*00 


0.0000D*00 
0.1500D*03 
0.0000D*00 
0.1000D*01 
0.0000D*00 
0.0000D*00 
-0.5000D*01 
0.0000D*00 
0.0000D*00 
0.0000D*00 
0.0000D*00 
0.0000D*00 
0.0000D*00 
0.0000D*00 
0.0000D*00 
0.0000D*00 
0.0000D*00 
0.0000D*00 
0.0000D*00 
0.0000D*00 
0.0000D*00 
0.0000D*00 
0.0000D*00 
0.0000D*00 
0.1703D*00 
0.0000D*00 
0.0000D*00 
-0.1703D*00 
0.0000D*00 
0.0000D*00 
0.1703D*00 
0.0000D*00 
0.0000D*00 
-0.17020%00 
0.0000D*00 
0.0000D*00 
0.0000D*00 
0.0000D*00 
0.0000D*00 
0.0000D*00 
0.0000D*00 
0.8035D*01 
0.0000D*00 
0.0000D*00 
0.0000D*00 
0.0000D*00 
0.0000D*00 
0.0000D*00 
0.0000D*00 
0.0000D*00 
-0.4033D*00 
0.0000D*00 
0.0000D*00 
0.6829D*00 
0.0000D*00 
0.0000D*00 
-0.1607D*00 
0.0000D*00 
0.00000%00 
0.15070“00 
0 0 
0.0000D*00 
0.1607D*00 
0.0000D*00 


0.0000D*00 
-0.1607D*00 


0.00000»00 
0.00000%00 
-0.18840+03 


0.00000+00 
0.0000D*00 
0.0000D*00 
0.0000D*00 
0.0000D*00 
0.0000D*00 
0.0000D*00 
0.0000D*00 
0.0000D*00 
0.0000D*00 
0.0000D*00 
0.0000D*00 
0.0000D*00 
0.0000D*00 
0.0000D*00 
0.0000D*00 
0.0000D*00 
0.0000D*00 
0.0000D*00 
0.0000D*00 
0.0000D*00 
0.0000D*00 
0.94200%02 
0.00000%00 
0.11150%01 
0.00000%00 
0.00000%00 
רת‎ 1 
0.0000D*00 
0.0000D*00 
0.1115D*01 
0.0000D*00 
0.0000D*00 
-0.1115D*01 
0 
0.00000%00 
0.00000%00 
-?- 523 
0.0000D*00 
0.0000D*00 
0.0000D*00 
0.0000D*00 
0.0000D*00 
0.0000D*00 
0.0000D*00 
0.0000D*00 
0.0000D*00 
0.0000D*00 
0.0000D*00 
0.00000+00 
0.00000%00 
0.00000%00 
0.00000%00 
0.00000%00 
0.00000%00 
0.00000%00 
0.00000%00 
0.0 
0.0000D*00 
0.0000D*00 
0.0000D*00 
0.0000D*00 
0.0000D*00 
0.0000D*00 


0.0000D*00 
-0.1500D*02 


0.0000D*00 
0.0000D*00 
0 .4 


-0.1500D*03 
0.0000D*00 
0.0000D*00 
0.135:D*01 
0.0000D*00 
0.0000D*00 

-0.6572D*01 
0.0000D*00 
0.0000D*00 
0.0000D*00 
0.0000D*00 
0.0000D*00 
0.0000D*00 
0.0000D*00 
0.0000D*00 
0.0000D*00 
0.0000D*00 
0.0000D*00 
0.0000D*00 
0.0000D*00 
0.0000D*00 
0.0000D*00 

- 12 
0.00000%00 
0.2306D*00 
0.0000D*00 
0.0000D*00 

-0.23040%00 
0.00000%00 
0.00000%00 
0.2306D*00 

-0.16300%00 
0.00000%00 

-0.2306D*00 
0.0000D*00 
0.0000D*00 
0.0000D*00 
0.64320%03 
0.00000+00 
0. 0000000 
0.00000%00 
0.00000%00 
0.00000%00 
0.00000%00 
0.00000%00 
0.00000%00 
0.00000%00 
0.0000D*00 
0.0000D*00 
0.0000D*00 
0.0000D*00 
0.0000D*00 
0.0000D*00 
0.0000D*00 
0.0000D*00 
0.0000D*00 
0.0000D*00 
0.0000D*00 
0.0000D*00 
0.0000D*00 
0.0000D*00 
0.0000D*00 

-0.1500D*02 
0.0000D*00 


0.0000D*00 
0.0000D*00 


0.0000D*00 
0.0000D*00 
- 6 


APPENDIX D 


COMPUTER PROGRAMS 
In order to plot the state variables of interest and compere them, the 
following computer programs were developed in WATFIV and used, in order 
to change the format of the OPTPLOT DATA, into a fixed format data file, 80 
characters in length with space between the values, to be used as input file 


for the EASYPLOT program. 


1. 0 Р e] ,.Op and ne 
$ JOB 
1 REAL A,B,C,D 
2 N= 1 
š WHILE (N.LE.100) DO 
4 READC5,10)A 
5 READC5,15)B 
6 READC5,20)C 
7 READC5,25)D 
8 WRITEC6,30)A,B,C,D 
9 10 FORMATCE1LS.7) 
10 15 FORMATC56X,E14.7) 
11 20 FORMATC42X,E14.7) 
12 25 FORMATCG2X,E14.7) 
15 50 ЕОБМАТ(8Ғ15.8) 
14 [+א=א‎ 
175 END WHILE 
16 STOP 
17 END 
$ENTRY 
2. Uncoupleg Roll Channel end for P, :? and $e 
$ JOB 
P СЕДЕА, ВСЮ, Е 
2 МЕ] 
5 WHILE CN.LE.100) DO 
4 05005, 70) 
5 READC5,15)B,C 
6 READC5,20)D 
7 READC5,25)E 
8 WRITEC6,30)A,B,C,E 
9 10 FORMATCE14.7) 
10 15  РОКМАТС13Х, 2Е1 4.7) 
11 20  ЕОРМАТСЕ] 4.7) 
12 25 FORMATCELG.7) 
1 50 ЕОКМАТСЧЕЈ 5.8) 
14 =א‎ 1 
1 5 END WHILE 
16 STOP 
17 END 
SENTRY 
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5. Coupled Pitch and Rol] Channels and for q, 2e. and nz 





$JOB 
1 REAL A, B, C, DBE- ‚б 
2 Ne] 
3 WHILE (N.LE.100) DO 
4 READ(5,10)A 
5 READC5,15)B 
6 READ(5,20)C 
7 READ(5,25)D 
8 READ(5,30)E 
9 READC5, 35)L 
10 READ(5,40)G 
11 WRITEC6,45)A,C,D,E 
12 10 ЕОВМАТСЕ14 7) 
1 5 15 FORMATCE14.7) 
14 20 FORMAT(56X,E14.7) 
15 25 FORMAT(G2X,E14.7) 
16 50 ҒОВМАТ(42Х,Е14.7) 
12; 55  FORMATCE14.7) 
18 40 FORMATCE14.7) 
19 45 ЕОВМАТСЧЕ15.8) 
20 N=N+1 
וק‎ END WHILE 
22 STOP 
23 END 
$ENTRY 
4 led P n | Channels and fer P, ? and 8e 
$JOB 
1 REAL A, B, C, D, E, F,G 
2 N=] 
3 HHILE (N.LE.100) DO 
ú READ(5,10)A 
5 READ(5,15)B 
6 READ(5,20)C 
7 READC5,25)D 
8 READ(5,30)E 
9 READ(5,35)F | 
10 READ(5,40)G | 
11 הזקון‎ О. 45)A,E,F,G | 
12 10 ЕОВМАТСЕ14.7) | 
15 15 FORMATCE14.7) | 
18 20 FORMATCE1G.7) | 
15 25 ЕОВМАТСЕ14.7) 
16 30 FORMAT(C56X, Elá.7) 
17 35 FORMAT(El4.7) 
18 40 ҒОВМАТ(42Х,Е14.7) 
19 45  FORMAT(4F15.8) 
20 N=N+l 
21 END WHILE 
22 STOP 
23 END 


SENTRY 
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APPENDIX E 


COUPLED STATE-FEEDBACK DESIGN INPUT DATA 
FOR POPLAR PROGRAM 
(UNIMPAIRED SYSTEM) 


1 1 
100.0 0.1 0,1 
1.0 0.00 0.0 
0.5 0.6 ol 0 
0 10 0 5 5 $ 1000 0 
22 "53 


-0.15000000000000%05 0.00000000000000%00 0.00000000000000%00 0.00000000000000%00 
0.0000000000000D*00 0.0000000000000D*00 0.0000000000000D*00 0.00000000000000%00 
0.0000000000000D*00 0.0000000000000D*00 0.0000000000000D*00 0.0000000000000D*00 
0.15000000000000%05 0.00000000000000%00 0.00000000000000%00 0.00000000000000%00 
0.00000000000000%00 0.00000000000000%00 0.00000000000000%00 0.00000000000000%00 
0.0000000000000D*00 0.00000000000000%00 0.0000000000000D*00 0.1553100000000D*01 
0.0000000000000D*00 0.1000000000000D*01 0.0000000000000D*00 0.0000000000000D*00 
0.0000000000000D*00 0.0000000000000D*00 0.0000000000000D*00 0.0000000000000D*00 
0.0000000000000D*00 0.0000000000000D*00 0.0000000000000D*00 0.0000000000000D*00 
0.0000000000000D*00 0.0000000000000D*00 0.0000000000000D*00 0.0000000000000D*00 
0.0000000000000D*00 0.0000000000000D*00 0.0000000000000D*00 0.0000000000000D*00 
0.0000000000000D*00 0.00000000000000%00-0.65720000000000%01 0.00000000000000%00 

-0.50000000000000%01 0.00000000000000%00 0.00000000000000%00 0.00000000000000%00 
0.0000000000000D*00 0.0000000000000D*00 0.00000000000000%00 0.0000000000000D*00 
0.0000000000000D*00 0.00000000^0000D*00 0.0000000000000D*00 0.0000000000000D*00 
0.0000000000000D*00 0.0000000000000D*00 0.0000000000000D*00 0.0000000000000D*00 
0.0000000000000D*00 0.0000000000000D*00 0.0000000000000D*00 0.0000000000000D*00 
0.0000000000000D*00 0.0000000000000D*00 0.0000000000000D*00 0.0000000000000D*00 
0.0000000000000D*00 0.1000000000000D*01 0.0000000000000D*00 0.0000000000000D*00 

-0.5926000000000D*02 0.0000000000000D*00 0.0000000000000D*00 0.0000000000000D*00 
0.0000000000000D*00 0.0000000000000D*00 0.0000000000000D*00 0.0000000000000D*00 
0.0000000000000D*00 0.0000000000000D*00 0.0000000000000D*00 0.0000000000000D*00 
0.0000000000000D*00 0.0000000000000D*00 0.0000000000000D*00 0.0000000000000D*00 
0.0000000000000D*00 0.0000000000000D*00 0.00000000000000%00-0.44360000000000%02 

-0.14920000000000%00 0.00000000000000%00 0.00000000000000%00 0.17887000000000%01 
0.00000000000000%00 0.00000000000000%00 0.00000000000000%00 0. 00000000000000%00 
0.00000000000000%00 0.00000000000000%00 0.00000000000000%00 0.00000000000000%00 
0.00000000000000%00 0.00000000000000%00 0.00000000000000%00 0.00000000000000%00 
0.0000000000000D*00 0.0000000000000D*00 0.00000000000000%00 0.00000000000000%00 
0.00000000000000%00 0.00000000000000%00 0.00000000000000%00 0.00000000000000%00 
0.0000000000000D*00 0.1000000000000D*01-0.7019000000000D*00 0.0000000000000D*00 
0.0000000000000D*00 0.0000000000000D*00 0.0000000000000D*00 0.0000000000000D*00 
0.0000000000000D*00 0.0000000000000D*00 0.0000000000000D*00 0.0000000000000D*00 
0.0000000000000D*00 0.0000000000000D*00 0.0000000000000D*00 0.0000000000000D*00 
0.00000000000000%00 0.0000000000000D*00 0.0000000000000D*00 0.0000000000000D*00 
0.00000000000000%00 0.00000000000000%00 0.00000000000000%00-0.44560000000000%02 

-0.14920000000000%00-0.10433600000000%04 0.00000000000000%00 0.00000000000000%00 
0.00000000000000%00 0.00000000000000%00 0.0000000000000D*00 0.00000000000000%00 
0.00000000000000%00 0.00000000000000%00 0.00000000000000%00 0.00000000000000%00 
0.00000000000000%00 0. 00000000000000%00 0.00000000000000%00 0.00000000000000%00 
0.00000000000000%00 0.00000000000000%00 0.00000000000000%00 0.00000000000000%00 
0.00000000000000%00 0.00000000000000%00 0.00000000000000%00 0.00000000000000%00 

-0.18840000000000%05-0.96200000000000%02 0.96200000000000%02-0.94200000000000%02 
0.94200000000000%02 0.00000000000000%00 0.00000000000000%00 0.00000000000000%00 
0.00000000000000%00 0.00000000000000%00 0.00000000000000%00 0.0000000000000D*00 
0.00000000000000%00 0.00000000000000%00 0. 00000000000000%00 0.00000000000000%00 
0.25042000000000%00 0.11154000000000%01 0.170520000000000%00-0.11154000000000%01 

-0.17020000000000%00 0.00000000000000%00 0.00000000000000%00 0.10092000000000%02 

-0.14300000000000%00 0.00000000000000%00 0.00000000000000%00 0.00000000000000%00 
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0.00000000000000+00 0.0000000000000D*00 
0.00000000000000+00 0 0 
0.00000000000000*+00- 00 0 
0.11154000000000*01 0 


0.0000000000000D*00 0.0000000000000D*00 
0.0000000000000D*00 0.0000000000000D*00 
0.0000000000000D*00 0.0000000000000D*00 
-0.1115400000000D*01-0.1703000000000D*00 


0.00000000000000+00 0.0000000000000D*00-0.1009200000000D*02 0.0000000000000D*00 


0.0000000000000D*00 0.0000000000000D*00 
0.0000000000000D*00 0.0000000000000D*00 
0.0000000000000D*00 0.0000000000000D*00 
0.2304200000000D*00 0.1115400000000D*01 


-0.14630000000000D*00 0.0000000000000D*00 
0.0000000000000D*00 0.0000000000000D*00 
0.0000000000000D*00 0.0000000000000D*00 
0.0000000000000D*00 0.0000000000000D*00 


0.17020000000000%00-0.11155000000000%01-0.17050000000000%00 0.0000000000000D*00 


0.0000000000000D*00 0.0000000000000D*00 
0.0000000000000D*00 0.0000000000000D*00 
0.0000000000000D*00 0.0000000000000D*00 
0.0000000000000D*00 0.0000000000000D*00 


0.0000000000000D*00 0.1009200000000D*02 
-0.1430000000000D*00 0.0000000000000D*00 
0.0000000000000D*00 0.0000000000000D*00 
0.0000000000000D*00 0.0000000000000D*00 


0.00000000000000%00-0.23042000000000%00-0.11154000000000%01-0.17030000000000%00 


0.00000000000000%00 0.00000000000000%00 


0.11154000000000%01 0.1703000000000D*00 


-0.1009200000000D*02 0.0000000000000D*00 0.0000000000000D*00 0.0000000000000D*00 
-0.1430000000000D*00 0.0000000000000D*00 0.0000000000000D*00 0.0000000000000D*00 
0.0000000000000D*00 0.0000000000000D*00 0.0000000000000D*00 0.0000000000000D*00 
0.0000000000000D*00 0.0000000000000D*00 0.0000000000000D*00 0.0000000000000D*00 
0.0000000000000D*00 0.0000000000000D*00 0.0000000000000D*00 0.0000000000000D*00 
0.0000000000000D*00 0.0000000000000D*00-0.1000000000000D*01 0.8871240000000D*03 
0.4432110000000D*03-0.4432110000000D*03 0.4432110000000D*03-0.44321100000000D*03 
0.0000000000000D*00 0.0000000000000D*00 0.0000000000000D*00 0.0000000000000D*00 
0.0000000000000D*00 0.0000000000000D*00 0.0000000000000D*00 0.0000000000000D*00 
0.0000000000000D*00 0.0000000000000D*00 0.0000000000000D*00 0.0000000000000D*00 
0.0000000000000D*00 0.0000000000000D*00 0.0000000000000D*00 0.0000000000000D*00 
0.0000000000000D*00 0.0000000000000D*00 0.0000000000000D*00 0.0000000000000D*00 
0.0000000000000D*00 0.0000000000000D*00 0.00000000060000D*00 0.0000000000000D*00 
0.00000000000000+00 0.00000000000000+00 0.00000000000000*:00 0 
0.0000000000000D*00 0.0000000000000D*00 0.0000000000000D*00 0.0000000000000D*00 


0.0000000000000D*00 0.0000000000000D*00 


0.0000000000000D*00 0.8034620000000D*01 


0.0000000000000D*00 0.0000000000000D*00 0.0000000000000D*00 0.0000000000000D*00 
0.0000000000000D*00 0.0000000000000D*00 0.0000000000000D*00 0.0000000000000D*00 


0.0000000000000D*00 0.1000000000000D*01 


0.0000000000000D*00 0.0000000000000D*00 


0.0000000000000D*00 0.0000000000000D*00 0.0000000000000D*00 0.0000000000000D*00 
0.0000000000000D*00 0.0000000000000D*00 0.0000000000000D*00 0.0000000000000D*00 
0.0000000000000D*00 0.0000000000000D*00 0.0000000000000D*00 0.0000000000000D*00 
0.0000000000000D*00 0.0000000000000D*00 0.0000000000000D*00 0.0000000000000D*00 
0.0000000000000D*00 0.0000000000000D*00 0.0000000000000D*00 0.0000000000000D*00 
0.0000000000000D*00 0.0000000000000D*00 0.0000000000000D*00-0.1760000000000D*02 


0.0000000000000D*00 0.0000000000000D*00 
0.0000000000000D*00 0.0000000000000D*00 


-0.8000000000000D*01 0.0000000000000D*00 
0.0000000000000D*00 0.0000000000000D*00 


0.0000000000000D*00 0.0000000000000D*00 0.0000000000000D*00 0.0000000000000D*00 
0.0000000000000D*00 0.0000000000000D*00 0.0000000000000D*00 0.0000000000000D*00 
0.0000000000000D*00 0.0000000000000D*00 0.0000000000000D*00 0.0000000000000D*00 
0.00000000000000+00-0.75325000000000+00-0.88352000000000*:00 0 


0.0000000000000D*00 0.0000000000000D*00 


-0.5000000000000D*01 0.0000000000000D*00 


0.00000000000000%00 0.00000000000000%00-0.40554000000000%00 0. 00000000000000+00 


0.0000000000000D*00 0.0000000000000D*00 
0.0000000000000D*00 0.0000000000000D*00 
0.0000000000000D*00 0.0000000000000D*00 
0.0000000000000D*00 0.0000000000000D*00 
0.0000000000000D*00 0.0000000000000D*00 
0.0000000000000D*00 0.0000000000000D*00 
0.0000000000000D*00 0.0000000000000D*00 
0.0000000000000D*00 0.0000000000000D*00 
0.0000000000000D*00-0.1027100000000D*00 
0.14518200000000%02-0.16181840000000%02 
0.0000000000000D*00 0.0000000000000D*00 
0.0000000000000D*00 0.0000000000000D*00 
0.0000000000000D*00 0.0000000000000D*00 
0.0000000000000D*00 0.0000000000000D*00 
0.10271000000000%00 0.12048000000000%00 
0.1418184000000D*02 0.0000000000000D*00 
0.0000000000000D*00 0.1406700000000D*00 
0.0000000000000D*00 0.0000000000000D*00 
0.0000000000000D*00 0.0000000000000D*00 
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0.0000000000000D*00 0.0000000000000D*00 
0.0000000000000D*00 0.0000000000000D*00 
0.0000000000000D*00 0.0000000000000D*00 
0.0000000000000D*00 0.0000000000000D*00 
-0.60000000000000%01 0.0000000000000D*00 
0.0000000000000D*00 0.6829100000000D*00 
0.0000000000000D*00 0.0000000000000D*00 
0.0000000000000D*00 0.0000000000000D*00 
0.0000000000000D*00 0.0000000000000D*00 
-0.120%8000000000%00 0.4795000000000D-01 
-0.1500000000000D*02 0.0000000000000D*00 
-0.1406700000000D*00 0.0000000000000D*00 
0.0000000000000D*00 0.0000000000000D*00 
0.0000000000000D*00 0.0000000000000D*00 
0.0000000000000D*00 0.0000000000000D*00 
-0.4795000000000D-01-0.1431820000000D*02 
-0.1500000000000D*02 0.0000000000000D*00 
0.0000000000000D*00 0.0000000000000D*00 
0.0000000000000D*00 0.0000000000000D*00 


0.0000000000000D*00 0.0000000000000D*00 0.0000000000000D*00 0.0000000000000D*00 
0.00000000000000%00 0.10271000000000%00 0.12048000000000%00-0. 47950000000000-01 
-0.14318200000000%02 0.14181840000000%02 0.00000000000000%00 0.00000000000000%00 


0.14067000000000%00 0.00000000000000%00 


0.0000000000000D*00 


-0.15000000000000%02 


0.0000000000000D*00 0.0000000000000D*00 0.0000000000000D*00 0.0000000000000D*00 
0.0000000000000D*00 0.0000000000000D*00 0.0000000000000D*00 0.0000000000000D*00 
0.0000000000000D*00 0.0000000000000D*00 0.0000000000000D*00 0.0000000000000D*00 
-0.10271000000000%00-0.12048000000000%00 


0.47950000000000-01 0.14318200000000%02 
0.0000000000000D*00 0.0000000000000D*00 


0.0000000000000D*00 


-0.1418184000000D*02 


-0.1500000000000D*02-0.1406700000000D*00 0.0000000000000D*00 0.0000000000000D*00 


0.0000000000000D*00 0.0000000000000D*00 
0.0000000000000D*00 0.0000000000000D*00 
0.0000000000000D*00 0.0000000000000D*00 
0.0000000000000D*00 0.0000000000000D*00 


0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 


0.0000000000000D*00 
0.0000000000000D*00 
0.00000000000000%00 
0.0000000000000D*00 


0.26988300000000%04-0.26988300000000%06-0.26988300000000%04 0.26988300000000%04 


0.00000000000000%00 0.00000000000000%00 
0.0000000000000D*00 0.0000000000000D*00 
0.0000000000000D*00 0.0000000000000D*00 


0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 


0.0000000000000D*00-0.1353100000000D*01 0.0000000000000D*00 


0.0000000000000D*00 0.0000000000000D*00 
0.0000000000000D*00 0.0000000000000D*00 
0.0000000000000D*00 0.0000000000000D*00 
0.0000000000000D*00 0.0000000000000D*00 
0.0000000000000D*00 0.0000000000000D*00 
0.0000000000000D*00 0.0000000000000D*00 
0.0000000000000D*00 0.0000000000000D*00 
0.0000000000000D*00 0.0000000000000D*00 
0.0000000000000D*00 0.0000000000000D*00 
0.0000000000000D*00 0.0000000000000D*00 
0.0000000000000D*00 0.0000000000000D*00 
0.0000000000000D*00 0.0000000000000D*00 
0.0000000000000D*00 0.0000000000000D*00 


0.0000000000000D*00 
0.00000000000000%00 
0.00000000000000%00 
0.00000000000000%00 
0.00000000000000%00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 


0.0000000000000D*00-0.9420000000000D*02 0.9420000000000D*02 


0.0000000000000D*00 0.0000000000000D*00 
0.0000000000000D*00 0.0000000000000D*00 
0.0000000000000D*00 0.0000000000000D*00 
0.0000000000000D*00 0.0000000000000D*00 
0.00000000000000%00-0.14300000000000%00 
0.23042000000000%00 0.00000000000000%00 
0.0000000000000D*00 0.0000000000000D*00 


0.94200000000000%02 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 


0.00000000000000%00-0.14300000000000%00 0.0000000000000D*00 


0.0000000000000D*00 0.0000000000000D*00 
0.00000000000000%00 0.0000000000000D*00 
0.23042000000000+00 0 
0.00000000000000%00 0.00000000000000%00 
0.44321100000000%03-0.46321100000000%03 
0.0000000000000D*00 0.0000000000000D*00 
0.0000000000000D*00 0.0000000000000D*00 
0.0000000000000D*00 0.0000000000000D*00 
0.0000000000000D*00 0.0000000000000D*00 
0.0000000000000D*00 0.0000000000000D*00 
0.0000000000000D*00 0.0000000000000D*00 
0.0000000000000D*00 0.0000000000000D*00 
0.0000000000000D*00 0.0000000000000D*00 
0.0000000000000D*00 0.1760000000000D* 02 
0.0000000000000D*00 0.0000000000000D*00 
0.0000000000000D*00 0.0000000000000D*00 
0.0000000000000D*00 0.0000000000000D*00 
0.0000000000000D*00 0.0000000000000D*00 
0.0000000000000D*00 0.0000000000000D*00 
0.0000000000000D*00 0.0000000000000D*00 
0.0000000000000D*00 0.0000000000000D*00 
0.0000000000000D*00 0.0000000000000D*00 
0.0000000000000D*00 0.0000000000000D*00 
0.00000000000000%00-0.12048000000000%00 
0.0000000000000D*00 0.0000000000000D*00 


0.0000000000000D*00 
0.0000000000000D*00 


-0.1884000000000D*03 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*0Q 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.6572000000000D*01 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 

-0.9420000000000D*02 
0.0000000000000D*00 

-0.14300000000000%00 

-0.23042000000000%00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 

-0.23042000000000%00 
0.00000000000000%00 


0.00000000000000%00-0.145300000000000%00 


0.00000000000000%00 


-0.44221100000000%03 


0.00000000000000%00 
0.00000000000000%00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.8835200000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 


0.0000000000000D*00 
0.4432110000000D*03 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 


0.0000000000000D*00-0.1500000000000D*02 


0.12048000000000%00 
0.00000000000000%00 
0.00000000000000%00 
0.0000000000000D*00 
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0.00000000000000* 00 
0.0000000000000D*00 
-0.1500000000000D*02 
0.0000000000000D*00 


0.0000000000000D*00 0.0000000000000D*00-0.1500000000000D*02 


0.0000000000000D*00 


0.0000000000000D*00-0.1206800000000D*00 0.0000000000000D*00 0.0000000000000D*00 


0.0000000000000D*00 0.0000000000000D*00 0.0000000000000D*00 
0.0000000000000D*00-0.1500000000000D*02 0.1206800000000D*00 
0.0000000000000D*00 0.0000000000000D*00 0.0000000000000D*00 
0.269883100000000+06-0.26988300000000+04-0 6 
0. . 0 Р 

0.0000000000000D*00 0.0000000000000D*00 0.0000000000000D*00 
0.0000000000000D*00 0.0000000000000D*00 0.0000000000000D*00 
0.0000000000000D*00 0.0000000000000D*00 0.0000000000000D*00 
0.0000000000000D*00 0.0000000000000D*00 0.0000000000000D*00 
0.0000000000000D*00 0.0000000000000D*00 0.0000000000000D*00 
0.0000000000000D*00 0.0000000000000D*00 0.0000000000000D*00 
0.0000000000000D*00 0.0000000000000D*00 0.0000000000000D*00 
0.0000000000000D*00 0.0000000000000D*00 0.0000000000000D*00 
0.0000000000000D*00 0.0000000000000D*00 0.0000000000000D*00 
0.0000000000000D*00 0.0000000000000D*00 0.0000000000000D*00 
0.0000000000000D*00 0.0000000000000D*00 0.0000000000000D*00 
0.0000000000000D*00 0.0000000000000D*00 0.0000000000000D*00 
0.0000000000000D*00 0.0000000000000D*00 0.0000000000000D*00 
0.0000000000000D*00 0.0000000000000D*00 0.0000000000000D*00 
0.0000000000000D*00 0.0000000000000D*00 0.0000000000000D*00 
0.0000000000000D*00 0.0000000000000D*00 0.0000000000000D*00 
0.0000000000000D*00 0.0000000000000D*00 0.0000000000000D*00 
0.0000000000000D*00 0.0000000000000D*00 0.0000000000000D*00 
0.0000000000000D*00 0.0000000000000D*00 0.0000000000000D*00 


0.0000000000000D*00 0.0000000000000D*00 0.0000000000000D*00 
0.0000000000000D*00 0.0000000000000D*00 0.0000000000000D*00 
0.0000000000000D*00 0.0000000000000D*00 0.0000000000000D*00 
0.0000000000000D*00 0.0000000000000D*00 0.0000000000000D*00 
0.0000000000000D*00 0.0000000000000D*00 0.0000000000000D*00 
0.0000000000000D*00 0.0000000000000D*00 0.0000000000000D*00 
0.0000000000000D*00 0.0000000000000D*00 0.0000000000000D*00 
0.00000000000000%00 0.0000000000000D*00 0.0000000000000D*00 
0.0000000000000D*00 0.0000000000000D*00 0.0000000000000D*00 
0.0000000000000D*00 0.00000000000000 “00 0.00000000000000“00 
0.0000000000000D*00 0.0000000000000D*00 0.0000000000000D*00 
0.0000000000000D*00 0.0000000000000D*00 0.0000000000000D*00 
0.0000000000000D*00 0.0000000000000D*00 0.0000000000000D*00 
0.0000000000000D*00 0.0000000000000D*00 0.0000000000000D*00 
0.0000000000000D*00 0.0000000000000D*00 0.0000000000000D*00 
0.0000000000000D*00 0.0000000000000D*00 0.0000000000000D*00 
0.0000000000000D*00 0.0000000000000D*00 0.0000000000000D*00 
0.0000000000000D*00 0.0000000000000D*00 0.0000000000000D*00 
0.0000000000000D*00 0.0000000000000D*00 0.0000000000000D*00 
0.0000000000000D*00 0.0000000000000D*00 0.0000000000000D*00 
0.0000000000000D*00 0.0000000000000D*00 0.0000000000000D*00 
0.0000000000000D*00 0.0000000000000D*00 0.0000000000000D*00 
0.0000000000000D*00 0.0000000000000D*00 0.0000000000000D*00 
0.0000000000000D*00 0.0000000000000D*00 0.0000000000000D*00 
0.0000000000000D*00 0.0000000000000D*00 0.0000000000000D*00 
0.00000000000000+00 0.0000000000000D*00 0.0000000000000D*00 
0.00000000000000+00 0.00000000000000+00 90 
0.0000000000000D*00 0.00000000000000%00 0.00000000000000“00 
0.0000000000000D*00 0.0000000000000D*00 0.00000000000000 “00 
0.0000000000000D*00 0.0000000000000D*00 0.0000000000000D*00 
0.0000000000000D*00 0.0000000000000D*00 0.0000000000000D*00 
0.0000000000000D*00 0.0000000000000D*00 0.0000000000000D*00 
0.0000000000000D*00 0.0000000000000D*00 0.0000000000000D*00 
0.0000000000000D*00 0.0000000000000D*00 0.0000000000000D*00 
0.0000000000000D*00 0.0000000000000D*00 0.0000000000000D*00 
0.0000000000000D*00 0.0000000000000D*00 0.0000000000000D*00 
0.0000000000000D*00 0.0000000000000D*00 0.0000000000000D*00 
0.0000000000000D*00 0.0000000000000D*00 0.0000000000000D*00 
0.0000000000000D*00 0.0000000000000D*00 0.0000000000000D*00 
0.0000000000000D*00 0.0000000000000D*00 0.0000000000000D*00 
0.0000000000000D*00 0.0000000000000D*00 0.0000000000000D*00 
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0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.26988300000000%04 
0.10000000000000“01 
0.00000000000000“00 
0.00000000000000“00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.1000000000000D*01 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.1000000000000D*01 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.1000000000000D*01 


0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.1000000000000D*01 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.1000000000000D*01 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.1000000000000D*01 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.1000000000000D*01 
0.0000000000000D*00 
0.0000000000000D*00 
0.00000000000000“00 
0.0000000000000D*00 
0.0000000000000D*00 
0.1000000000000D*01 

0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 

4.00000000000000%00 
0.0000000000000D*00 

0.1000000000000D*01 

0.0000000000000D*00 
0.0000000000000D*00 

0.0000000000000D*00 
0.0000000000000D*00 

0.0000000000000D*00 


0.00000000000000%00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 


0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.00000000000000%00 
0.00000000000000%00 
0.00000000000000%00 
0.00000000000000%00 
0.00000000000000%00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 


0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.00000000000000%00 
0.00000000000000%00 
0.00000000000000%00 
0.00000000000000%00 
0.00000000000000%00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.00000000000000%00 
0.00000000000000%00 
0.00000000000000%00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.00000000000000%00 
0.00000000000000%00 
0.00000000000000%00 
0.00000000000000%00 
0.00000000000000%00 
0.00000000000000+00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
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0.1000000000000D*01 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.10000000000000%01 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.1000000000000D*01 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.1000000000000D*01 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.1000000000000D*01 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.1000000000000D*01 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.1000000000000D*01 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.1000000000000D*01 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.1000000000000D*01 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.1000000000000D*01 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.1000000000000D*01 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.0000000000000D*00 
0.1000000000000D*01 
0.0000000000000D*00 


পারা... чаи 
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